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CALCULATION OF POWER SPECTRUMS FROM DIGITAL TIME 


SERIES WITH MISSING DATA POINTS 

INTRODUCTION 

Quite frequently vyhen analyzing time series, gaps are encountered in the data. For example, 
time records of three-day averages of the brightness temperatures of the earth’s surface made by 

■ pv . 

a microwaii j radiometer onboard an eart -monitoring satellite contain a number of missing data 
arrays (Wilheit 1972; Zwally and Gloersen 1977; Rayner and Howarth 1979). In the analysis of 
tides, observations occur randomly in time (Zetler 1965). In astronomy and oceanography, factors 
such as weather and the availability of measuring instruments dictate when observations can be 
made. Even with portions of a data record missing, there may be a random element in the spacing 
of the points. Unfortunately, standard formulas for calculating power spectrums (in order to 
detect periodicities in the data) are based upon the measurements or observations being equally 
spaced. 

The problem of estimating the power spectrum for unequally spaced data has been studied by 
various authors. A number of articles and papers dealing with the subject have appeared in the 
literature within the last ten years. 

Meisel ( 1 978, 1979) discusses the problem of unequal spacing and describes procedures which 
interpolate in the time domain. He also describes a general procedure for processing segmented 
data sets and for obtaining the Fourier transform of a time series with arbitrary spacing. His pro- 
cedure is most efficient for data where large gaps occur between the data sets and the data sets 
are made up of equal or nearly equal intervals (within 20%) of time. An excellent set of references 
on the subject of unequally spaced data can be found in both of his articles. 

Deeming (1975) studied the Discrete Fourier Transform of a function defined for arbitrary 
data spacing and also discussed properties of spectral windows for various kinds of data spacing. 


Vanicek (1969, 1971) combined least squares fitting procedures and Fourier analysis in 
order to remove imdesired inniiences on the spectrum, such as biases, linear trends, etc. His method 
can also be used to uncover periodicities from unequally spaced time series. Other references 
which discuss the least squares aspect of the problem are Baming (1963), Taylor and Hamilton 
(1972), Lomb (1976), Ferraz-Mello (1977), and Wells and Vanicek (1978). 

Ferraz-Mello (1980) defines a “Date Compensated Discrete Fourier Transform” for estimating 
the power spectrum. His estimate is obtained by orthonormalizing the three basis functions — 
cos cot, sin cot, and 1, using a Gram-Schmidt procedure. Since the set of basis functions is com- 
plete, harmonic filtering of the time series is made possible. 

The present analysis treats the problem of missing data from the autocorrelation function 
(or lagged products) point of view. It is assumed that there are gaps in the data. However, data 
points which are separated by gaps are assumed to be nearly equally spaced (to within approxi- 
mately 20% of an equal sampling interval). 

Two algorithms are developed for calculating the power spectrum from the autocorrelation 
function when there are missing data points. Both procedures use an average sampling interval 
in order to compute the lagged products where the sampling interval has been calculated over sets 
of data points separated by gaps. The intervals between all adjacent data points are then calcu- 
lated to the nearest multiple of this average spacing. Two counters are set up. One counter sums 
the number of lags of a given size, the other the corresponding cross product terms. 

One procedure, the Correlation Function Power Spectrum, takes the Fourier transform of 
the lagged products function directly to obtain the spectrum, and is equivalent to the Fourier 
Transform Power Spectrum method when there are only missing data points (no random or arbi- 
trary spacing of points). The other procedure, the Modified Blackman-Tukey Power Spectrum, 
takes the Fourier transform of the mean lagged products (obtained from the ratio of the two 
counters mentioned above) to calculate the power spectrum. 
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Both algorithms are compared with the Fourier Transform Power Spectrum method (Deeming 
1975) and two least squares procedures (one by Vanicek 1971 and the other by Ferraz-Mello 1980). 

Examples are given showing recovery of frequency components by all the techniques from 
simulated periodic data where portions of data are missing and where random noise has been added 
to both the time points and to the values of the time series. In addition, the methods are compared 
using real data. All procedures performed j^ually well in detecting periodicities in the data. 

Power Spectrum From the Autocorrelation Functions 

Let the time series be given by the pairs of points (tj, for i = 1,2,. . .N. 

The discrete Fourier transform Xf,^(a)) of x^ is defined by (Deeming 1975 and Appendix B) 

N 

i=l 

= (0l-j02)x (1) 

where 


0 -( 01 , 02 ) = 



and 


x"^ = (Xi,Xj, . . .,Xj„) 


The normalized Fourier Transform Power Spectrum Sp(cj) is then (Appendix C) 

^ Nx^x 

and the Least Squares Power Spectrum is (Vanicek 1971 and Appendix C) 
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( 3 ) 
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In the case of equal spacing the expressions in (2) and (3) are equivalent since 

(0^0) = (N/2) I with I the 2 X 2 

identity matrix 


The real symmetric matrix <fxpJ can be written as 

00"^ = I + C + 

where I is the (N X N) identity matrix and 



Since 

x‘*'(C + C'^)x = 2x'^Cx 
the spectrum in (2) can also be written as 

2 

Sp(co) = NxTx + 2x^Cx) 


(4) 


(5) 


( 6 ) 
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or as 


Now Jet 


Sp (co) 


2 

NxT'x 



x-Xj coscoAtg - 



yk = XiXj i=l,2, ...,N 

Atj^ = At^ j = 1,2, . . „N 

k= 1,2,. . •.N(N+l)/2 


( 8 ) 


(9) 
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and define y_|^ and At_j^ by 


y-k = ^j^i = yk 

At_k = Atji 

cos(cjAtj^) = cos(cjAt_jj) 


( 10 ) 


Then the expression in (8) can be written as the discrete Fourier transform of the set of cross 
product terms yj^ 




2 

Nx^ X 
2 

Nx'^x 


lc=-K 

K 

2 cosCwAt^) 


( 11 ) 


k= -K 


where K = N(N+l)/2. 


Let At be an average spacing calculated over data points which are separated by gaps and let 
Atjj be an approximation to Atj^ obtained by rounding Atj^ to the nearest multiple of At. 

At. 


At^ 

Ar. = !L + 0.5 At 

k L 


(12) 


and [x] is the largest interger contained in x. 


Then (1 1) can be approximated by the Correlation Function Power Spectrum S(.(co) defined 


as 


M 


i - S ‘'k cos(kojAt) - xT) 


(13) 


k=0 


where the lagged product terms Cj^ are given by 


K 




(14) 


i=0 


•5(Ar.) = 


1 for At- = kAt 


I 0 for Atjt^ kAt 


(15) 


and 


and M, the maximum lag in the data is 


M = 



( 16 ) 


If there are only gaps in the data and the points are separated by multiples of the average 
spacing. At, then (13) will be the same as the Fourier Transform Power Spectrum in (2). Further- 
more, if the maximum lag M is replaced by m where m is 70% to 80% of M, the resulting spectrum 
will be about the same as (13). 


It should be noted that for equally spaced data, using the Blackman-Tukey procedure (Black- 
man 1958) mean lagged products Rj. are first calculated 


N-r 

1=1 




(17) 


(i=0,l, . . . ,m) 


Then the spectrum is given by 




(18) 


We can modify this standard procedure of Blackman and Tukey by defining two counters 




i=0 

- 5 ] yi5(Arj) 


(19) 


i=0 


where 6(Ar-) has been previously defined in (15) above and y^ in (9) 


Then, mean lagged products c^ will be given by 

c'='k 

and we will have a modified Blackman-Tukey spectrum Sjyjg j(cj) 


( 20 ) 
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( 21 ) 


I 

= -Y, Ck™5<klJ^t)--xTx 

V '^"0 

Comparing (17) and (21) we see that the counters nj^ and ensure that when the time differ- 
ences across all time points are calculated (multiples of the average time difference or spacing), the 
cross products are summed with their proper lag number k. No interpolation is done and only 
known data points are used. Also, no zeros are substituted for missing data points. Thus, for a lag 
of k, instead of N-k points in the sample there will be something less. If not too many points are 
missing, the c|^ terms will not be greatly affected. For a value of m equal to about 60% of the maxi- 
mum lag M, the results obtained by using (21) will be about the same as ( 13). It should be noted 
that in computing j(co) some negative amplitudes will occur which is due to the fact that we 
are dividing by a number less than or equal to N-k (instead of N) and we are not using all the lagged 
product terms. 

Frequency Aliasing 

If a period function, cos(27rft), is sampled at equally spaced intervals of time. At, the following 
. values of the function are obtained 

Xj = cos(27rfiAt) (22) 

(i=...-2,-l,0,l,2,. ..) 

A sampling interval of At seconds gives rise to a maximum resolvable (cutoff) frequency of 
fj. = (l/2At) Hz. The function Xj can be written as 

Xj = cos|27T[(2kf^, - OiAt]} (23) 

(k=l,2,. ..) 

(i=... -1,0,1...) 

From (23) it can be seen that for |2kf^ - f | < f^., the frequencies |2kf^, - f| and f are indis- 
tinguishable and are called aliases of each other. 


EXAMPLES 


In order to demonstrate the algorithms developed in (13) and (21), data points Xj 
erated using the following function 

X. = sin(27rfi tj) + cos(27rf2 tj) + sin(2:rf3 tj) 

+ cos(2trf4 tj) + sin( 2 jrfs tj) 

where fj = 0.37 Hz, fj = 0.14 Hz, fj = 0.42 Hz, f 4 = 0.65 Hz, and fj = 1.09 Hz. 


were gen- 


( 22 ) 


First 101 points were calculated at equal intervals of 1 second with tj = i — 1 (i=l,2,. . .101). 
With this choice of sampling interval aliasing can be expected. This was done in order to assess the 
effects of the sample interval upon recovery of the frequency components. 


The power spectrum was then calculated by the Correlation Function method in (13) with M 
(the maximum lag in the data) replaced by m = 0.6M, the Modified Blackman-Tukey procedure in 
(21) with M replaced by m = 0.5M, the Fourier Transfer method (Deeming 1975) in (2), and the 
Least Squares procedure (Vanicek 1971) in (3). The results can be seen in Figures 1 through 4. 
Note that the frequencies 0.09 Hz and 0.65 Hz which show up in the spectrums are aliases of 1.09 
Hz and 0.65 Hz. These particular frequencies, since they are greater than the Nyquist or cutoff 
frequency (f^ = 0.5 Hz), beat with twice the cutoff frequency. Thus, (fs-2fc) and (2fc-f4) show up 
in the spectrum. 


Next, approximately 30% of the data points were deleted. The spectrum for each method was 
calculated and can be seen in Figures 5 through 8. Note that the frequency components have been 
recovered (except for aliasing). 

Figures 9, 10, 11, and 12 demonstrate recovery of frequency components when random noise 
from a Gaussian distribution with zero mean and standard deviation of 0. 10 has been added to the 
time points tj. The function in (22) was calculated for these pseudo-random time points and the 
spectrum obtained using all four methods. It should be noted that with a standard deviation of 
0. 1 (and a sampling interval of ) second), approximately 95% of the noise will be less than or equal 
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to 0.2 or two sigmas. This is in agreement with Meizel (1978) who concluded that as long as data 
spacing was within 20% of equal intervals ordinary sine [(sin x)/x] interpolation is satisfactory. 

Finally, Figures 13, 14, 15, and 16 show recovery when there are both missing data (the same 
points being deleted as in Figures 5 through 8) and random noise (the same noise as in Figures 9 
through 12). 

In Appendix D there are additional examples using the same time series but with a sampling 
interval of 0.72 seconds, noise having a standard deviation of 0.2 seconds on the time points and 
amplitude noise with a standard deviation of 20% of the value of the function. Also the method of 
Ferraz-Mello (1980) is compared with the above methods. All of the procedures produced com- 
parable results and the frequencies in the data were recovered equally well. 

As a further example to demonstrate both the Correlation Function Power Spectrum method 
in (13) and the ^Modified Blackman-Tukey spectrum in (21) and compare them with other spectrum 
calculations, we have used real data. 

A covariance matrix was built up from 1 70 three-day average brightness temperature maps 
(Wilheit 1972; Zwally and Gloersen 1977;Rayner and Howarth 1979) of the earth’s surface where 
the temperatures were obtained from an Electrically Scanning Microwave Radiometer (ESMR) on- 
board an orbiting NIMBUS 5 satellite. The raw data was averaged in both time and space. The 
time period extended from September 1973 through May 1975 and the spatial data consisted 
of an area in the South Polar region approximately 3780 kilometers by 4380 kilometers. Land 
masses were masked out. Thus the observations are essentially the brightness temperatures of 
the water and sea ice. Out of a total of 201 possible three-day average maps covering the time 
period, there were 3 1 missing data arrays. An eigenvalue/eigenvector analysis of the covariance 
matrix was made and the principal components calculated. The first ten principal components 
accounted for more than 91% of the total variance in the data. These principal components (vary- 
ing in time across the 170 time'points) were then normalized. The second principal component 
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time series is plotted in Figure 17. Figures 18 through 25 are normalized power spectrums which 
have been calculated by various methods and which show the same periodicities in the data. The 
first peak at a frequency of about 0.0028 cycles/day corresponds with a period of about 360 days. 
The second peak of about 0.006 cycles/day corresponds with a period of about 167 days while 
the third peak at 0.0083 cycles/day corresponds with a period of about 120 days. These particu- 
lar periods can be related to the time variation of the spatial mean latitude of the outer ice boun- 
dary (360 days), the warm season (l67 days) and the spring summer removal of ice (120 days) 
reported by Rayner and Howarth (1979). Figure 18 is the power spectrum computed by least 
squares using (3). Figure 1 9 is the power spectrum computed by the Fourier Transform method 
using (2). Figure 20 is the spectrum computed by the method of Ferraz-Mello (1980). Figure 21 
has been calculated using (13) and is exactly the same as the Fourier Transform spectrum as there 
are only gaps in the data. Figure 22 has been calculated using instead of M, a value of m equal to 
80% of M and using (13). Figures 23 and 24 are for values of m equal to 70% and 60% of M respec- 
tively. Finally, Figure 25 has been calculated by the Modified Blackman-Tukey Power Spectrum 
algorithm in (21) with the lag number m equal to 60% of M. 

Analysis of Figures 18 through 25 show that all of the methods perform about the same. They 
differ only in the amount of calculation ne.cessary to obtain the spectrum. Using the Correlation 
Method in (13) where the lagged products are summed to about 70% of the maximum lag M and the 
Modified Blackman-Tukey procedure in (21) with an m equal to about 60% of the maximum lag M 
involves less computation and computer time than the other procedures. 

SUMMARY 

Two algorithms have been developed for calculating power spectrums from the autocorrela- 
tion function for time series with missing data points. Both methods use an average sampling 
interval to compute the lagged products. One procedure, the Correlation Function Power Spec- 
trum, calculates the Fourier transform of the lagged products directly to obtain the spectrum. This 
method is equivalent to the Discrete Fourier Transform spectrum (£)eeming 1975) when there are 
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only missing data points (that is, the time series is equally spaced except for gaps) and all lagged 
products are included in the computation. The other procedure, the Modified Blackman-Tukey 
Spectrum, calculates the Fourier transform of the mean lagged products in order to obtain the 
spectrum. Satisfactory results can be obtained by summing lagged product terms up to approxi- 
mately 50% to 80% of the maximum lag in the data. 

The algorithms have been compared with the Discrete Fourier Transform method (Deemin ; 
1975) and two least squares procedures (Vanicek 1971 and Ferraz-Mello 1980) using both simu-' 
lated and real data. In the case of the simulated data, there are portions of the time series missing 
and random noise with a Gaussian distribution has been added to both the time points and to values 
of the function. All of the methods performed satisfactorily and frequency components in the data 
were recovered equally weU. 

Recovery of frequency components in the presence of noise having a standard deviation of 
10% of the sampling interval (Gaussian noise with 95% less than 20% or 2 standard deviations) 
supports Meisel (1978) that as along as the data spacing is within 20% of equal intervals, ordinary 
sine [(sinx)/x] interpolation is satisfactory. 
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APPENDIX A 


Definitions 

A real-valued function x(t) is periodic with period T if for all t 

x(t + T) = x(t) 

The Fourier transform X(f) of x(t) is defined as 


X(0 



x(t)e“j^’’^^* dt 


A.1 


A.2 


if the integral exists for all f, where f is frequency and t is time . 

If x(t) is periodic and limited over some interval T^, then x(t) can be represented by its Fourier 
Series 


where f^ = 1 /T^ and 


x(t)= 2^ ttjje 

n=-~ 


j.27rn fgt 


.To/2 

M “J ^ « ii V 

“n = I , x(t)e ° dt 


I 


-j 27 rn f-t 


-To/2 


A.3 


A.4 


It should be noted that in taking the Fourier transform of a periodic function x(t), the Fourier 
coefficients which are obtained, are the same as the Fourier transform evaluated at n/T^. 
Another useful relationship is the correlation integral given by 


oo 

z(t) = J x(t) h(t + T)dr 


A. 5 


It can be shown (Brigham 1974) that the Fourier transform Z(f) is given by 

Z(0 = H(f) -X^Cf) A.6 

where X*(f) is the complex conjugate of X(0- 
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The discrete analogue of A. 5 and A.6 is that if 

N-i 


z(k) = x(i) • h(k + i) 
i=0 


A. 7 


(k = 0, 1,2,...N-1) 

is the correlation function for x(i) and h(i), then the Fourier transform Z(n) of z(k) is given by 

Z(n) = X*(n) -H(n) A. 8 

where the discrete Fourier transform X(n) of x(k) is defined as 


N-l 

X(n)= x(k)e-j 2 ’T"k/N 
k=0 

(n = 0,l,...,N-l) 

If x(t) and h(t) are the same function in A.5, we have the autocorrelation function of x(t) 


A.9 


z(t) = 1 x(r)x(t + T)dr 

or if x(i) is the same as h(i) in A.7, the discrete autocorrelation function of x(i) is 

N-l 

z(k) = y x(i)x(i + k) 
i=0 


A.10 


A.11 


(k = 0. 1,...,N-1) 

Then, using A.6 and A. 8, the Fourier transform of the autocorrelation function Z(f) is given by 

Z(0 = X*(f) •X(f) = | X(f)(2 A.12 

in the continuous case and by 

Z(n) = X*(n) -X(n) = | X(n)|2 A.13 

in the discrete case. 

Equations A.12 and A.13 represent the contribution of the frequency f (or n) to the total power 
in the waveform. 
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Tlie total power in the time domain is equivalent to the total power in the frequency domain 
by Parsevalls equation which is given by 


OO 



x^(t)dt = 



A. 14 


in the continuous case and by 


N-l N-1 

2^x2(k) = (I/N)2^X2(n) 

k=0 n=0 


A.15 


in the discrete case. 

The relationship in A.l 2 (or equivalently A. 13) is an important one in that it provides us with 
another tool for obtaining the power spectrum of a time series. Note that the power spectrum 
can also be obtained by taking the square of the magnitude of X(0 in A. 2 or of X(n) in A.9. 
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APPENDIX B 


Definition of Impulse Function and Properties 

The impulse function 5(t) is a distribution which assigns to a test function x(t) the value of 

x(0) 


J 


5(t)x(t)dt = x(0) 


The function 5(t - t-) is defined as 


< 

f 


5(t-to)x(t)dt = x(to) 


From the above we obtain the following 


f 


8(t)dt= 1 

The product o f a 5(t) fu nction by an ordinary function y(t) is defined by 

oo 

[5(t)y(t)]x(t)dt = J 6(t)[y(t)x(t)] dt 




and if x(t) is continuous at t = ri 


B.l 


B.2 


B.3 


B.4 


5(to)x(t) = 5(to)x(to) 


B.5 


Derivation of Discrete Fourier Transform for an arbitrary spacing 

Let x(t) be the function to be sampled and define the sampling function s(t) as a series of im- 
pulse functions to be applied at the times tj^ 


( 1 ) 


The symbol 5(t) for the impulse function is not to be confused with the symbol 5 used in the text. The 5 there 
is the Kronecker delta function. 



Then the sampled function is 


s(t) 


IN * 

= ^5(t-t^) 


k = l 


N 

x(t) = x(t)^^ 
k=l 


The discrete Fourier transform of x(t) is then 


Xn(0 = 



x(t)e"^^^^^ dt 



5(t - t]^)e“j^’^^‘ dt 


Interchanging the integral and summation in B.8 and using B.2 and B.5 

N 

k=l 


r 6(t-tj^)x(t)e-j2Jrftjjt 
^ _oo 


=z: 


x(t^)e 


-jZTTftk 
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APPENDIX C 


MATHEMATICAL DEVELOPMENT OF 
LEAST SQUARES POWER SPECTRUM 

The standard. least squares problem of determining a best fit of the data to an assumed mathe- 
matical model is equivalent to a minimum norm problem in Hilbert space. The latter approach is 
taken by Vanicek ( 1 969, 1971) and Wells and Vanicek ( 1 978) in their treatment of the least squares 
power spectrum. 

In this appendix we will develop the optimum least squares estimate of the power spectrum 
from both viewpoints. A geometric interpretation of the problem can be seen in Figure 1. 

Figure 1 

Hilbert space H 

Let X be any (measurement) vector in a Hibert space H where x is given by 

x'^ =(Xi,X2, . . . , Xjj) C.l 

and where Xj is the value of the time series at the time tj. 

Let M be a closed subspace spanned by the (Nxl) vectors and <^2 

( coscuti sincoti 

cosu?t2 sino;t2 

costotf,, sincotpi 

Then, the projection theorem assures us that there is a unique N vector p = 0$ such that^ ^ ^ 

Ilx-Pll < ||x-p|| for all p e M C.3 

where the norm ||x|| is defined in terms of the inner product of x with itself 

l|x||^=x^x C.4 

^The use of the caret in p and c refers to the least squares estimate and is not to be confused with Wells and Vani- 
cek’s (1978) use of this symbol which is used to refer to the space $ of “constituents” (bias term, linear trend, 
quadratic trend, etc.). We are*assuming no constituents and therefore Vanicek’s p = 0. 
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The necessary and sufficient condition that p Exists is that (x-p) be orthogonal to M, or, in 
particular 0. That is, the inner product of f x-p) and <j> must be equal to zero. 


0 ^(x- 0 c) = 0 


C.5 


Equation (C.5) is the vector representation of the usual normal equations. Note that since 
and 02 span the space M, they are linearly independent, and therefore, the Gram matrix 0 0 is non- 
singular. Hence, we can solve for c uniquely. 

C = (0^0)~*0^X C .6 


Then, in the least squares sense, the best approximating vector to x (which is contained in M) 
is the projection of x onto the subspace M. See Figure 1. Therefore, 

p = 0c C.7 

This particular vector gives the minimum norm ||x-pll or the distance from x to the subspace M. 


An optimum normalized least squares spectrum 8 ^( 0 ;) can now be defined [Wells and Vanicek 
(1978)1 as 


Sl(w) = 


llx-pll^ 

llxll* 


C.8 


A “nice” feature of this particular choice for the spectrum is that it has a minimum of zero 
when the vector x is orthogonal to M (or when the projection of x onto M, p, is equal to zero), and 
a maximum when x is contained in M (that is, when the projection of x onto M is equal to x and 
therefore (x-p) = 0). 


Using the fact that 


p^(x-p) = 0 


and the definition of the norm in (4) we can write 


T'' 
X p 

Sl(^) ~ X 
^ x^x 


C.9 
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CIO 


or using C.7 


^ ^ ^ X^0((/>^0) '0^X 

S, (cj) Y 

^ X* X 


C.ll 


In the usual least squares terminology we wish to fit data points (tj.Xj) to some assumed 
mathematical model 

p(t) = cos(o)t) a + sin(a;t)b C. 1 2 

determining the least squares estimates a and B for a and b such that 

N 


Q= ^ |xj- [(coso)tj)a + (sincotpb]| 


C.13 


i=l 


is minimized. This is equivalent to finding p to minimize ||x-p|| in C.3. 

The standard procedure is to differentiate C.13 with respect to both a and b, then set these 
expressions equal to zero. Two equations in two unknowns result 


^ cos* ojtj 1 ^ ^ ^ (sinwtjXcoscutj) | b = ^ Xjcoscutj 

i=i / \ 

N \ ^ \ ^ 

(coscjtjXsincJtj) I a+| sin*wtj Jb = Xjsincjtj 


C.14 


i=l 


i=l 


i=l 


which are the normal equations or in matrix notation equation C.5. 


The two elements of the vector £ in C.6 are the estimates t and B and we have 

c'^ = (a,B) Cl 5 

We therefore obtain, as above, the best estimate for p (or p) in M as equation C.7. 

It is worthwhile at this point to compare the expression obtained for the least squares spec- 
trum with the Discrete Fourier Transform spectrum, Sp(a)). 

Recall that the Discrete Fourier Transform Xf,j(cj) of Xj (Appendix B) is given by 
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X^(u,)= Y, 

i=l 

N N 

= ^ Xjcos(cjtj) - j ^ XjSin(o)tj) 
i=l i=i 

In matrix notation since 0i and are defined in C.2 and x in C. 1 

Xj^(O») = 0i' x-j02X 

= (0T-j^I)x 


C.17 


Then the amplitude or power spectrum will be given by 

11X^(00)11^ =X;li(co)TxN(cu) 

= x'^(0, +j02> (0T“J>T)X 

= X^00^X 

where X|j(oj) is the complex conjugate of Xj^j(oo). 


We can now write a normalized Fourier Transform Spectrum as 


^ . ,_2x'^00'‘'x 
Sp(w) — 

Nx‘ X 


C.18 


C.19 


where this latter expression is obtained by noting that in the case of equally spaced data we will 
have 

(^T^ = (N/2)I C.20 

where I is the 2X2 identity matrix. 


In the case where the data is equally spaced the least squares estimate S^fco) in C. 1 1 will be 
equivalent to C. 19, the Discrete Fourier Transform power spectrum. 


Finally, it should be noted that if the spacing is not too irregular (within 20% of equal 
sampling) the two spectrums will be about the same. 
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MODIFIED BTACKMAN-TUKSY AND CORRETATTON FUNCTION PO’'fER SPECTpU?.'! ^POC-PAV 
FOR UNSQUAT TY SPACED- DATA POINTS 


00610 C 

000c;0 

00010:0 

000401 

000601 

00070 C 
0001:501 C 
0001901 C 
001001 
001101 
001 £01 
001 1301 
0101401 
001501 
00 1601 
001701 
001:30 
001901 
002010 
00210 
00££0i 

010230 

00240 

00250 C 

00£60i C 

00£70 

00£:30i 

00£90i 

003001 

00310 C 

0013E01 

003301 

0013401 

0035:1" 

00 360 
0013701 
003:30 
00390 
00400 
0014101 
00i4£0i 
01014301 

0vi44i3 

00450 
0046'0 
004701 C 
0014:301 C 
004901 C 
0050115 C 
00510 0. 
0i05£0i C 
005301 C 
0015“'" C 
00550 C 


MRIH riF:Ii..iER FOR FOICGR 
IMPLICIT RERL»3 (Pi-HjO-Zj 
D I MEHS I ON X ( 5010 ) 5 T ( 500 )- j RH 1 ( 50i0i ) 

DIMENSION NT I M 1 : 500:1 
D I MENS I ON TM ( 500i :• j IKM ( 5Gi0i :i > RN ( 50i0i :i 
DIMENSION DEU:500:i 

:5ETTING UP SIMULPlTED DRTfl (T ( I :i j K t I j :i , 1 = 1, NT I'iHERE 
Ti:i:i CONTRINS RRNDOM NOISE RNi;i;i, ZERO MERN RND STRNDRRD 
DEiJIRTION OF 0. IG FROM R GRUSSIRN DISTRIBUTION 
NT = 139 

FI = 3. 14159£67D0i 
RRGl = £3 0iD0*FT '^0L 37D0 
RRG£ = £3 0DG-«FT*0L 14D0 
RRG3 = £ ^ 0D0i^P I *01 * 4£D0i 
RRG4 = £.0D0i*PI*03 65D0 
RRG5 = £ ^ 0D0i*P I * 1 * 019D01 
DO 1 I = 1,NT 
T 1 : 1 :i = 01 . 7£D0* ( I - 1 :i 

RN 1 : 1 :i = ERRN 1 1 : - 1 , 01 , 1 £7:37 , 0 - 0D0 ,0-1 0D0 :i 
Ti:i;i = Ti:i:i + rni:i:i 

1 : 1 :i = DS I N ( RRG 1 *T ( I ) :• •• DCOS ( RRG£*T ( I :i :i + DS I N 1 : RRG3*T ( I :i :i + 
iDCOS i:rrg4*t ( I :i :i +dsin i:rrg5*t i: i ;i ;i 
STDEU = 0- 10D0*Xt:i) 

RN 1 1 : 1 ;i = BRRN 1 1 : - 1 , 01 , 1 37 : 37 , 0 - 0 D 01 , STDE'.' :i 

RDDING GRUSSIRN NOISE MITH ZERO MERN RND SIGMR EQURL TO £0i7. OF 
THE RMPLITUDE OF THE FUNCTION- 
X(i;i = X(i:i + RNKD 
1 CONTINUE 
DO £ I = 1,NT 
£ NTIM 1 : 1:1 = 1 

PUTTING 'HOLES’ IN THE DRTR 
DO 4 I = 1,8 
NTIM 1:1 + 10:1 = 0 
NT I MCI + 40:1 = 0 
NTIM Cl +70 j = 0 
4 NTIM C 1 + 90:1 = 0 
L = 0 

DO 3 I = 1,NT 
IFCNTIMCi:i -EQ- 1:1 L = L + 1 
iFCNTiMci;i ,EQ- 1:1 tmcl:i = tci:i 
3 IFCNTIMCD -EQ- 1:1 XMCU = XCIj 
NT = L 

CRLL PONCOR C 0 , 0D0 , NT , 0i - 60D0 , 0 - 000 1 DO , 1 00i , XM ' TM 'i 

:3TOP 

END 

SUBROUTINE PONCOR CSi-J, NT, FRLRG, FST, Ni-J, F, T:i 

SI'i IS R :3NITCH PRRRMETER- IF SN = O MODIFIED BLRCKMRN-TUKEV 
SPECTRUM IS CRLCULRTED- IF SN = 1, THE CORRELRTION FUNCTION 
POICR SPECTRUM IS CRLCULRTED - 
■ NT IS THE TOTRL NUMBER OF POINTS 

TCD RRE THE TIME POINTS RND FCn RRE THE i.CRLUES OF THE TIME 
:3ERIES Cl = 1, NT) - 

M IS THE MRXIMUM LRG USED FOR COMPUTRTION PURPOSES RND IS 
EQURL TO R FRRCTION FRLRG OF THE MRXIMUM LRG IN THE DRTR MRXLRC 
NN IS THE NUMBER OF FREQUENCIES CRLCUuRTED EETICEN THE STRRTINC 




FREQUEhCV FST RND THE CUTOFR FREQUENCY 
iNPLii-IT ^.ERiL*b 


UUS'rU 

C0R.00 


:jT ppTr:TC pi Tv M|j^■ T T.-i! Jf .1 | hC ThDEY 

li i PlENU i Uh h I. btiC j ? T ( b00 J ? C ■; 500 j 
D I MENS I ON SM 1 1 5oti ) .•> SME ( 500 j j TM ( 500 ) j FM ( 500i J 


F16.9) 


I,- 


Fijh LL-r1R‘! ! i H ! j ■ ;rTR!_ r'L^RF’L*'T 


001610 C 

OBTAINING AUERRGE TIME 

0101680 

CALL DELTAT ( T j NT j 4 j 1 . 0 

0101630 

MRITE (4!. 999) DELT 

0101640 

999 FORMAT ( IX d DELT = , 

0101650 

MMl = M - 1 

0101660 

PI =- 3. 141598653590D01 

001670 C 

SETTING COUNTER EQUAL 

01016801 

DO 3 1 = 1? 5010 

006901 

SMld) = 0 

01017001 

SM£d) = 0 

0017101 

3 C ( I ) = 01 

0017801 

NMl = NT - 1 

0107301 C 

REFERENCING TIME POINT- 

010740 

TMd) = 0.0D0 

0107501 

DO 4 I = 1?NM1 

0107601 

DELA = Td + 1)-Td) 


00770 4 tmi;i+i:i = tm(i;i + delh 

010 7801 DO 5 I = 1 j NT 

00790 '5 T i; I :i • = TM i I j 

00S00 C ROUrffilNG TIME POINTS TO THE NERREST MULTIPLE OF R.'.JERfiGE TIME 

00S10 C DIFFERENCE DELT OBTAINED FROM DELTRT SUBROUTINE 

00:i£0 DO 6 I = IjNT 

00830 MTM = TM ( I ) ••••■DELT + 0.5D0 

00840 6 TMi;i:i = MTM 

00850 C ELIMINATING DUPLICATE TIME POINTS. 


010860 

L = 0 

008701 

DO 7 I = 1?NM1 

008801 

IFCTMd) .EQ.TM'd+D) GO 

008'?0 

L = L + 1 

010'?001 

Ti;L) = TMd) 

00-? 10 

FN!:L) = F(I) 

00'?80 

7 CONTINUE 


T(L+1) = TM(NT) 

010940 

FMi'L+D = F(NT) 

0109501 

NT = L 

0i0i‘560 

NMl = NT - 1 

0i0'r‘70 

MMl = M-1 

0i0i‘980i C 

CALCULATING THE MAXIMUM 


00990 C 
0110100 ( 
01010 
01080 
010301 
Vlll040 
011060 
01060 
01070 
0101801 
010901 C 
011010 C 
01110 C 
01180 C 
01 130 
011401 
01150 
01160 
01170 


USE M = FRLAG:«MA:KLAGj LAGS OF A GI'.'EN SIZE AND CORRESPONDING 
SUNS OF LAGGED PRODUCTS. 

MAKLAG = T (NT) - T(l) + 1.0D0 
M = FRLRG^-^MAXLAG 
DO 10 I = 1 j NT 
DO 9 J = I, NT 
I LAG = T (J) - Td) + 1.0D0 
SM 1 ( I LAG ) = SM 1 ( I LAG ) + 1 . 0D0 
9 SM8(ILAG) = SMSdLAG) + FMd)«FM(j;i 
10 CONTINUE ■ 

CALCULATING LAGGED PRODUCTS OR MEAN LAGGED PRODUCTS FROM 
COUNTERS SMI (I) RND SM£d) AABO'.'E. 

IFlSI'J.EQ. 1.0D0) GO TO 18 
DO 11 I = IjM 

I F ( SM Id) . LE J 1 . 0iD0 ) C ( I ) = V0 . 0D01 
I F ( SM Id) . LE . 1 . 0iD0 ) GO TO 11 
C(I) = SM8d).-SMl d) 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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01 180 
01 180 
01800 
01810 
01880 
01830 
01840 
01850 
01860 
01870 
61880 
01880 
013010 
01310 
01380 
01330 
01340 
013501 
01360 
01370 
01380 
013801 
014010 
01410 
01480 
01430 
01440 
01450 C 
01460 C 
01470 C 
0148:0 
01490 
015001 
01510 
01580 
01530 
01540 
01550 
011560 
01570 
015:30 
015801 
01600 
01610 
01 1 6801 
016301 
01640 
01650 
01660 
01670 
0116:501 
01680 
0170101 
01710 
01780 
01730 
017401 
01750 
017601 
01770 
017:50 
0 1 780 


11 COMTIhUE . • 

GO TO 14 ' 

18 DO 13 I = IjM 

13 0(1) = Sr'18(I) 

14 FHORM = 0G0II0 

DO 15 I = 1j NT 

15 FNORM = FNORM +. FMi:i;i^«^«8 

FEND = 1.0D0i.--i;8.0D0^«DELT:i 
DELF = (FEND - FST) .••■DFLORlT i;NI-J;i 
FREQ = FST 

DO 18 I = 1.- NN 
SUMR = 01.0D0 
DO 16 J = l 5 N 
FJ = J - 1 

RRG = 8.0D0»PI>FR.EQ»:F.J«DELT 

16 SUMR = SUMR + C C J) :«DCOS (RiRG) 

IF CSI'UEQ. 1.0D0) GO TO 17 

SPOOR = 800.0D0«(8.0D0»:SUMR - FNORM.-'DFLOfiT (NT) ) .•■ FNORM 
GO TO IS 

1 7 SPOOR = 800 . 0D0* ( 8 . 0D0^SUMR-FNORM ) ( DFLORT ( NT ) *FNORM ) 
1:5 I'iRITE (45 80 J FREGis iPOOF: 

18 FREQ = FREQ + DELF 
80 F0RMRT(8F16.8) 

8888 RETURN 
END 

RERL FUNOTION BRRN1«8 ( 1 5 IKCT 5 IFRN 5 RMERN 5 SD) 

I MPL 1 0 1 T RERL^:5 ( R— H 5 Oi~Z 5 :f ) 

:5D THE DESIRED STRNDRRD DEU 

RNERN-..THE DESIRED MEAN. 

H. ..THE POP SIZE. 

DRTR -IHERE.■••■187S7.•••• 

DflTP^H.•■•36. 0D01.-' 

IF (I KEY) 55454 

4 I HERE = IFRN 

5 IF i;i) 65757 

6 0:RLL GRUSS ( IHERE 5 :5D5RMEfiN5'..iRL5H;t 
IFRN = I HERE 

GO TO 8 

7 0:RLL RRNDU C IHEF:E5 IFRN 5 '..'flU 
IHERE = IFRN 

8 BRRNl = I..IRL 
RETURN 

END 

:5UB:R0UTINE GRUSSaX5S5RM5'..'5H:i 
IMPLIOIT REALMS (R-H5 0-Z5 :t) 

K = H 
R = 01.0D0 
DO 50 I = I 5 K 
OREL RRNDU(IX5 IY 5 V) 

IX = lY 
50 R = R + Y 
HO = H.--18.0D0 
H8 = H.--8.0D0 

I. .' = (S:«(R-H8) ).-DSQRT(HO;i + RM 
RETURN 

END 

SUBROUT I NE RRNDU ( I X 5 I Y 5 YFL ) 

I MPL 1 0 1 T RERL:«S ( R-H 5 0-Z 5 :f ) 

DATA .J.J.J5.--1087.-- 
lY = IX«.JJJ5 
IFdY) 55656 

5 lY = lY + 814748.3647+1 
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01 00-0 
W 1 1-! 1 0 


■0 1 CC0 


'01S30 


01 840 


01850 


01860 


01870 

1“ 

01880 

j- 

018”30 

c 

01900 

c 

01910 

1" 

01980 


0 1 930 

1“ 

01940 

c 


YFL = lY . . 

YFL = YFL», 4656b 13D-Y 

RETURh 

EfiD 

SUBROUTiriE l!ELTfiTi;T.N,KLOOPjTOLjDELT:i 
IMPLICIT REfiL»3 m-HjO-Zj 

n I MEr Hs I OM T i: 500 ) j del ( 500 ) » belm 1 500 j 

TCI) = TIME POINTS FROM I = 1j H 

KLOOP = fi PfiRfiMETER BEFINING THE NUMBER OF TIMES THE 

SUBROUTINE GOES THROUGH H LOOP TO THROH OUT TIME BIFFERENCE 
GREATER THAN TOL^^SB 

SB = STANBARB BEiJIATION OF THE TIME BIFFERENCES 

BELi\ = A‘.iERAGE TIME BIFFERENCE CALCULATEB FROM SaJBFOJTINE 

XM i Ri.iERAGE TIME BIFFERENCE. BELT TAKES THIS 

iJALiJE UPON E:aT FROM THE SUBROUTINE. 


0195£i 


NMl = 

N - 

1 

01960 


KL = 1 



0 1 970 


BO 1 I 

— 

IjNMI 

011980 

1 

BEL ( I f 

— 

T(i+i;i - T(i;i 

0 1 9901 


SUMl = 

0 


1.130001 


SUM8 = 

0 


08010 


BO 8 1 

= 

1 J riMl 

018080 


SUMl = 

:5L 

Ml + BEL (I) 

08030 


:i;!JM8 = 

SL 

M8 + BEL(i:i»»8 

080401 


m = 3 

UMl 

.•■ BFLOAT (NMl ;i 

0i80i50i 


RRG = 

( SIJM8.-BFL0AT (NMl ) J - 

080601 


IF (RRG 

.LE 

. 0 . 0B0 j SB = 0 

080701 


IF (SB. 

EQ. 

0.0B0) GO TO 999 

08011501 


SB = B 

SORT (ARGJ 

080901 


IF (SB. 

LE. 

YM) GO TO 999 

081001 


L = 0 



081101 


BO 3 I 

= 

1j NMl 

081801 


IF (BEL 

(i;i 

.le.tol«sb;i l = l 

08130 


IF (BEL 

(D 

.le.tou«sb;i belm(l 

08140 


IF(KL 

.GT 

. KLOOP) GO TO 999 

08150 


KL = K 

L + 

1 

08160 


NMl = 

L 


08170 


BO 4 I 

— 

1 jNMl 

08180 

4 

BEL ( I i 

= 

BELM ( I j 

018190 


GO TO 

39 


088001 

9S9 

BELT = 

XM 


08810 


RETURN 



08880 


ENB 



ENB OF 

l-i 





+ 1 


= BELtn 


INPUT 

0ii:c30 

QEB 


COMPUTATION OF 

THE MOB I FI EB B 

FROM THE pond: 

R PRUURRM. Nul 

PLUS HOLES IN 

THE BRTR. 60i:-; 

BELT = 0.7 

1 6986ys3 

01. U001 10.00010 

-1.548176786 

0.007073818 

-0. 150789885 

0.014046484 

-0.569139551 

y 4 yiil 1 1 ot' 

-0.754189357 

01 . 087998848 

-0 . 07 1 476469 

01 . 0349660601 

-1.853636884 


lk sFElTkun I. bi'i — y.i 
RNB THE FUtiL-T I ijfi 
: USEB FOR M. 


^.38 


ti. 041939£?1 

-01. 70101 :' 016 :-'S 

fi . 

? , 1 


P 

y . 06E35S90? 

“i j 

0.0fc'rci3c:il9 

1 . bo333;" 6ot' 

01. 076S05331 

3 . 6367 1 3433 

01. 6S377S543 

-1 . 047745357 

01. 01901751755 

1 . 10176301133 

01 . 0197734967 

—01. 139533379 

0. 10469S179 . 

0 . 683379 1 73 

0. 111671390 

0.409943943 

0. 113644603 

1 . 333373473 

0. 135617314 

0.911733360 

0. 133591036 

1 . 396405503 

0. 139564333 

17.457179640 

01. 1465374501 

7 . 933633 :^69 

0. 153510663 

-3.471330306 

0. 160433374 

4.315563713 

01. 1674570136' 

—3 . 633339737 

0. 1744301393 

-0 . 36 1 30130136 

0. 1314013510 

-0. 179330473 

01. 133376731 

-1 . 333336937 

0. 195349933 

4.793633391 

0.303333145 

-3.333365853 

0.309396357 

3.377395993 

0. 316369569 

-1.384619373 

U 4 dcLScL^d. 'i' 1 

-3.339691359 

0. 330315993 

3. 043041039 

0.337139305 

-3.041714553 

0.344163417 

1 . 196875771 

0-351135639 

0. 149344333 

0.353103341 


01 . 365033053 

<1* 

0 - 373055364 

—5 . 4363378014 

0. 379033476 

4. 633737647 

01. 3360011638 

—8 . 669465699 

0.393974900 

1 . 0330196453 

0.399943113 

3Si . 9656596331 

0. 30693 1034 

13. 639i;t3 73001 

0.313394536 

-4.313917966 

01. O'iz.0lb6r' :’43 

-1 . 68895;- 33 1 

0 . 3373409601 

—3 . 366650455 

0.334314171 

3.333904671 

0.341737333 

-5. 167401777 

0 . 3437601595 

5. 5903015537 

0.355^3:';fi09 

-1 . 9;3:-;4,450S;-; 

0 . 363701701 19 

5 . 09904 1451 

01. 3696301331 

•46' . 5350060014 

0 . 376653443 

33. 143117579 

01 . 333636655 

-10.739409463 

0 . 3901599367 

-3.304033313 

0 . 3975730179 

7 . 3353701909 

0.404546391 

—3 . 7696501657 

0.411519503 

0.141104103 

0.413493714 

•j’o .i y yfcfc'-.' .•*' 



0. 433439133 

-11. 040359636 

01. 4o!941c:-350i 

6 . 898448010101 

0. 446335563 

—1 . 713196398 

0 . 453353774 

1 . 449988304 

0.460331936 

5. :-'45934093 

01 . 467305 1 93 

-3. 0143039330 


39 



Hi.4f-4£7®4ie 

5. 96-0431300 

0.431351633 

1 - S3SS3S335- 

0. 4SS;3346'33 

3.379054135 

0.495193045 

1 . 36^755033 

0. 503171357 

-4.439314707 

0.509144469 

01. 3801715653 

0.516117681 

-1 . 331533319 

01 . 5330i90iS93 

0 . 8693013336 

0 . 5300164 1 05 

-1 . 165751987 

0.53703731? 

-0i.S07'337399 

0. 5440il053’9 

3.407178831 

01. 5501983741 

-3. 5376163301 

0.557956953 

0.759 171 349 

01. 564930164 

-0 . 309 1 90605 

■01. 5719013376 

-3. 168750319 

0.57SS765S8 

1.791153931 

01. 5S5849S0101 

-1 . 768993475 

01. 593833013 

-01 . 01 1 901019763 

0 . 599796334 

-0. 18018501333 

01 . 606769436 

0 . 485933476 

01. 613743648 

01 . 9793956801 

01 . 630715860 

—3 . 0i786‘4'34b3 

01. 6376890173 

1.981811540 

01. 634663383 

-3. 485895741 

0 . 64 1 6''S5495 

0.998141407 

01 . 648603701? 

13. 1650157663 

01.655581919 

' 3.804976381 

0.663555131 

-3.087106153 

01. 669538343 

01. 689841 10i5 

0.676501555 

1.855003671 

0.68'3474767 

3.603868543 


e.6'?i344?97'? 1 - 6?3795c::53 

EHIi OF DRTR 


COMPUTRTIOH OF THE CORRELRTICN FUNCTION POI'JER :i 
FROM THE POI'JCOR PROGRRM. NOISE ON TIME POINTS 
PLUS HOLES IN THE DRTR. 70>. OF THE MRKIMUM LRG 


BELT = 0.71> 

:’936883 

01 . 00101 1 00010101 

—01 . 3455601677 

01. 01017073313 

01 . 5401355097 

0. 01 4046434 

—01 . 363687738 

0.031019636 

0.695394888 

0 . 037993848 

—0 . 4ic! 1 6773S-I6 

01 . 0134966060 

0.483970455 

0.041939371 

-0.338465377 

0.048913483 

3.536149713 

0.055885695 

1.551455083 

0 . 01638589017 

0. 1 13976701 

0.0698331 19 

1 . 173969943 

01. 0768015331 

1 . 553435450 

0 . 0183778543 

01 . S'3339 1 8301 

0 . 09075 1 755 

1.417395157 

01. 097734967 

1 . 01401930636 

01. 104698179 

01 . 797709866 

0. 111671390 

1. 173196338 

01. 118644603 

01. 584598515 

0. 135617814 

3.300831346 

0. 1335910136 

1 . 737363367 


;PECTRUM (Si’J = i;i 
RND THE FUNCTION 
USED FOR M. 


40 



0 . i 39564 c!3S 

1 01 . 17 6’666633' 

0- 144F;3^4=;m 


0- 153510633 

1 . 435407500 

0i. 160 i433!374 

1 , 464490436 

0!-'16?45?0S'6 

-0.334136615 

0. l?4430c9S 

01. 531557907 

0. 131403510 

-0 . 30101573396 

0. 13S376731 

1 . 766309932 

0. 135349933 

1 . 567139366 

y ^ cLydoci'j 1 "45 

■ 0.637719004 

01. 3019296357 

1 . 764335963 

0. 316369569 

0.363335954 

y - uiii: Jc!4c y 1 

0.460436377 

0.330315993 

3.553330647 

0. 3371S9305 

1 . 557343741 

0.344163417 

0. 593495341 

0.351135639 

3.436166153 

0.353 103341 

0.537155337 

0 . 3650330152 

0.596353564 

0 . 373055364 

01 . 0i40i3330i99 

0.379033476 

1.011 154596 

0 . 33600 1 633 

0.676407792 

0 . 393974900 

3.337354015 

0.399943113 

13.409046133 

0.306921334 

3.730617639. 

0.313394536 

01.315599416 

01. 330367743 

3.373113740 

0 . 337340960 

-01 . 554626704 

0.334314171 

3. 494545994 

0.341737333 

-1.129107705 

01 . 343760595 

3. 165733654 

0 . 355733307 

3.731719033 

0.3627070119 

5.935293015 

0.369630331 

34 . 775661369 

01. 376653443 

1 1 . 739464133 

01. 333636655 

-0.304^1 .301 13 

0 . 390599367 

1 . 20137'60666 

0 . 397573079 

3.335475343 

■0^ 40454g,c;9i 

3 . 370993375 

0.411519503 

1 . 5015307399 

0.416493714 

20. 133406537 

01 . 43546’5936 

9.017391033 

0.433439133 

-01.510713677 

0.439413350 

3. 145133533 

0 . 446335563 

-0. 143131663 

01. 453-353774 

3 . 5563 1 09 1 6 

0.460331936 

3. 364341999 

01 . 467305 1 93 

3.357335519 

0.474373410 

3. 167749053 

0.431351633 

c! 4 

0.433234333 

3. 137511333 

0.495193045 

01. 163113331 

0.503171357 

01. 163394933 

0.509144469 

0.370633053 

0.516117631 

0 4 7'354'?'4'?7'5 

01. 5330901393 

0. 5150113359 

0.530064105 

-0 . 006443373 

0.537037317 

0. 923734335 

0.544010539 

1.045559133 

0.550933741 

0. 137313995 

0. 557956953 

0. 734370916 

0.564930164 

0. 513339336 


-41 


0 - 5 1 9030 : ■■ 6 

0. 136901111 


0. 33961 4 94:i 


0.538331349 


01 . 01 01 6 = 94 c! 01 6 > 

0 ^ ^i997'3bii!ii.'4 

1 . l4o!3c:' 7c5 

0 . 600769436 

0.715353603 

0.613743648 

0. 363155934 

01. 6301715860 

0.361367736 

0.637689073 

01 . 669~0i7536 

01 . 634663383 . 

01 . 5450168484 

01.641635495 

1. 148651415 

01 . 6486018701? 

6.31 383016 1 6 

0.655581919 

3. 149508091 

0.663555131 

. 01.377665375 

0.6695383^ ; 

01. 9v*i 1879757 

0.67650155 !< 

1.385836394 

01. 683474767 

3. 0801647764 

01 . 690447979 

1 . 759819507 


EHD OF DFlTfi ‘ 
INPUT 


■'A V' 



JEAST SQUARES POWEP SPECTRUM PRO'^-PAM FOP 
UMEQUATTY SPACED CATA POINTS (VAN'^CEK) 


60010 IMPLICIT REALMS (fl-HjO-Z) 

000£0 C POWER: SPECTRUM BY LEAST SQUARES METHOD (i;ANICEK;i 
00030 C HT = HUMBER OF TIME POINTS 

V00040 C T(I) = TIME POINTS. ..I = 1? NT 

60050 r: :ki:i;i = i.ialues of time series at points t(i:i...i = i- nt 

000S0 C FST = STARTING i.iALUE OF FREQUENCY 
00070 C FEND = ENDING 'JALUE OF FREQUENCY 

000 R 0 r: FNI'J = NI'J = TOTAL NUMBER OF FREQUENCY POINTS IN SPECTRUM 

00090 C DELF = FREQUENCY INCREMENT 

00100 C FREQ = CALCULATED FREQUENCY 

00110 C SPLSQ = CALCULATED LEAST SQUARES POhiER SPECTRUM 
00130 DIMENSION T(500 )j X(500'.i 

00130 PI = 3.1415S367D6 


00140 


READ (3j1) NT 

60150 


FuRMAT 1 . 1 10.! 

00160 


DO 3 I = 1j NT 

00170 

c 

READ t :5 ? :5 ) T ( 1 .1 > X ( I ) 

00180 

:5 

F ORMAT ( 3F 1 6 . 9 ) 

00190 


XNORM = 0 

00300 


DO 4 I = li. NT 

00310 

4 

XNORM = XNORM + Ki:i)»*«3 

Vj0330 


READ 1:3 .- 5:1 FSTj FENDi. FNI-J 

00330 

5 

FORMAT 1 : 3D 15. 8 J 

00340 


DELF = (FEND - FST) .--FNIu 

00350 


FREQ = FST 

00360 


NW = FNI'J 

00370 


i 'JR I TE ( 4 j 6 ) FST > FEND » FNW » DELF 

003S0 

6 

FORMAT (5D15. 8) 

00390 


DO 91=1- NW 

00300 


SI =0 . .i 

00310 


S3 = 0 

'00:33v3 


S:3 = 0 

y0U30 


S4 = 0 

00340 


S5 = 0 

00350 


DO 7 J = 1j NT 

00:060 


ARG = 3.0D0«PI»FREQ:«Ti;J) 

00:570 


CC = DCOS ( ARG ) 



'•5S = DS I N 1 , ARG .1 

00:590 


SI = SI + 

004001 


:53 = :53 + i.:u^SS 

00410 


:5:5 = 3:5 + SS:*:*3 

00430 


84 = 84 + CC'^X(.J) 

004:50 

"7 

S5 = S5 + SS»^Xi.J) 

0044lTi . 


BET = 31*S3 - :S3*«3 

00450 


A = (S3*S4 - S3*S5) ••■•'DET 

00460 


B = i-S3»S4 + S1«S5)--DET 

00470 


:SPLSQ = 100. 0D0* (S4-^A + S5*B:) 

00480 


I'iRITE (4 j8) FR'.EQ’ :SPLSQ 

00490 

8 

F URMAT 1 . izFiiti . 0 .1 

00500 


FREQ = FREQ + DELF 

00510 

999 

STOP 

00530 


END 

END OF 

DATA 



;norm 


:norm 


43 


iT FROM THE LEhST $C 


'5 1 D-UE 6 , o049T05yD 

00 0. 1010000 

£K &Gt495Ci5C' 

y4S‘:7:'yy5 

01 , 0ti99titi9y 

0.03651585 

0.0 1485 149 

01. 44401031 1 

0. 01980198 

1.31714898 

0. 0£'4?5848 


0. 0£9?0£9? 

1 . 63098383 

Ci Ci'-:4^‘^":4^ 

iJH 1 Cl 1 

0 . 03960396 

1 . 80859705 

0. 04455446 

01. 83'168'46'9 

0.04950495 

1 . 88491303 

01. 015445545 

8.11 973686 

01. 05940594 

0. 501440354 

01. 06435644 

01. 058379101 

0 . 069301693 

0.34417895 

0. 017485743 

1.17334088 

y .i y 'rciy 

8. 35680947 

0.08415848 

S- ‘?6'?c963?* 

0. 089101891 

16. 0i50i86410i 

0 . 0940594 1 

8.96340748 

0 . 0990109901 

8. 30814413 

0. 103960401 

1 . 08946886 

01. 10891089 

0.65786654 

0.11386139 

01. 89481986 

0.11881188 

0. 59331941 

0. 18376838 

0. 63078881 

0. 18871887 

1 . 03833314 

01. 13366337 

6 . 46884938 

0. 13861386 

15. 04579567 

01. 14356436 

101. 38980558 

0. 14851485 

8.86773144 

0. 15346535 

1 . 05348949 

0. 15841584 

01 . 8383366 

01. 16336634 

0.01088108 

0.16831683 

1 0.33904558 

0. 17386733 

8. v380930i53 

0. 17881788 

8. 18344881 

0. 18316838 

y . 46848848 

0. 18811881 

8. 943301418 

0. 19306931 

0.77934555 

0. 19801980 

1.83183637 

01. 80897vli30 

W. •.'4:' ?8l1iI. 

01. £0179801^9 

0 . 336870701 

0.81887189 

0. 10467851 

y -cl : ' C=C‘ 1 C 

a I'l -;fl 1 Sfi 

y - CCC: : C!CC’ 

0 . 05 

y. 88778877 

6. 039458 

y - C --Ct* C-C 

0. 15911101 

y- cy :-y 

0.80££8£5£ 

0.84857486 

01. 1 018868 IH 

0.84758475 

0.87366385 

0.85847585 

>•_! - c ^ rc* * cb :■ 

0.85748575 

01. 0i6.;-i63444. 

01. 86837 684 

0.08951885 

0 . 86 :-'386 4 

01. 160’98£37 

■ y - C c'c' .-■* C C’ 

01. 875360189 

0.87788773 

8. 39888580 

y - C'l-C! 1 CCC 

1 . 49959306 

y - CC= 1 c! 

4.08516080 


.44 







0.8980:- 981 

■ 0.6:-860815. 


£i . 89701897 1 



0.301 98080 

8. 3S6S4938 


0.30693070 

. S 6 1 ‘ £i401SI£i 


0.31188119 

1.10876838 


0.31683169 

4.33145590 


0.38178818 

S.8016OiOi63 


0i. 38673868 

0. 49719433 


0. 33168317 

1.53448408 


0 . 33663367 ' 

0.40988499- 


0.34158416 

8. 78419688 


0 . 34653466 

17. 80938933 


0.35 1485 15 

81.17850843 

£ 

0 . 35643565 

4. 891601669 


0.36138614 

0 . 93567389 


Oi . ;-;6A3366.4 

14. 9756796'9 


0.37188713 

84. 18459074 



6. 71830568 


0.38118818 

5.60381791 


0i. 38613868 

3,686891014 


0.39108911 

1.11 46 1 963 

■ 

0 . 39601396 1 

0 . 8397 4468 


0. 400990 l£i 

4 . 06886947 


0. 405940601 

8.79730388 


0.41089109 

0. 14878981 


0.41584159 

6.74S83350 


0. 48079808 

88 . 066 1 9043 


£1 . 485 4c'5s 

8. 58463337 


0. 43069307 

0.88476587 


0. 43564357 

5.08176888 

, - - ' 

0. 44059406 

0.89688580 


0. 44554456 

0.51 ^36868 

‘ v- / 

0 . 45049505 

1 . 37887334 


0.45544555 

3.05885336 


0 . 460396014 

1.84985137 

'Sifr-'--' 

0. 46534654 

0 . 86564688 


0. 47089703 

01.096981017 


0.47584753 

: 0 . 84054058 


01.48019808 

0. 18096984 


0.48514858 

1 . 88183370 


0 . 49009901 1 

0.71376860 


0.49504951 

• 0. 80393343 


0 . 5OiOi0OivIiOi0 

0.73493889 



tHI' UF DFiTR 


INPUT FOR LEPiST S:QUflRES POiCR SPECTRUM i:<.iPlNICEK;i fiND DISCRETE FOURIER 
TRfiNSFORM PQI.JER SPECTRUM I'JMERE THE SPfiCIMC IS 1 SECOND EXCEPT FOR 

HOLES IN THE DfiTR 


£l 0 


C. VZ100000000 


1 > 00001000001 1 , 79£. 1 87857 

c ^ 00100001000 - 1 , 4339859 1 6 

3 . OiOi0OiOiOi0Oi0 8 ^ 7083 1 588 1 


4 , OiOiOiOiOiOiOi0Oi 

5 . Oi00Oi0Oi0OiOi 
fc, , 0OiOiOiOi000Oi 
7 . OiOi0OiOi0OiOi0 


-1.747774168 
-01.881831906 
1.953107-48 
-1 . 591861893 



45 



C» ‘IT* '-jI* CO CO CC* CC» CC* O'J CC' CC> CC» 00 '••J '••J tT* ’O'* O'* O'* 0« O'* O'* O'* O'* O'* CH CH CR C.H »in CO CH CH CH CH 4* 4^ O'. O'J OO O'J O'J C*j C*j C*j C*j C*j PC' PO PO PO P*.' P'j PO PC' PC' PC' 


S- 000000000 

0. 577521509 

9- 000000000 

-0- 510703131 

1 ^ 0000n00MM 

i QO-j -;1 “pi 

19 -000000000 

-2- 030147064 

20 . 0000000100 

1 . 53:-!5:-:0ifi-7£ 

2 1 . 0100101001010101 

-2. 192373373 

22 . 0010101001010101 

2.210496S9S 

23 . 0101000VZ1010101 

0. 6570996121 

24 . 0010101010101010 

-0. 304905312 

25 . 0101001001010101 

0.999999122 

26 . 0101010100100101 

-0. 150379346 

2? - 010100100100101 

0 . 44364 1 0501 

23 . 0101001001010101 

0 . 332475263 

29 . 0100101010101010 

0. 792349873 

30 - 000000000 

-1 . 642040024 

3 1 . 00101010101010101 

-0.603910431 

3ii: - 000000000 

-1.343346922 

33 . 010100010101001 

-1.669335104 

34 , 01010010100'0101 

1 . 740466054 

35 , 00101010101010101 

0.357962310 

36 * 01010101010100101 

2.7469663U7 

37 - 000000000 

1 - 01746744201 

3S - 000VI100000 

-0. 133607533 

33 ^ 0010010001101 

—01 . 5'0iS33i5466' 

43 < 0101001010101001 

y - cl c! t'C!‘C‘C‘c! 

49 ^ 0100101010101001 

2.003250345 

501 < 010010101000101 

-01 . 0100001 1 707 

5 1 ^ 010101010101010101 

0.4421647001 

52 . 0101001001010101 

-0.629492118 

53 ^ 0010010101010101 

-2. 4538759301 


^ 000000000 -1^921 432014 1 

, 001001010010101 01 < ??S?7 1 ?€.4 


56 . 001001001010101 

*“1 - lc!iilc!yc*c*'4 

57 . 000010101010101 

2.339344455 

53 . 010101Vll010100101 

^ 1 4-*4c'c!y5 

59 . 01010101000010 

-1 . 579394735 

601 . 0100101001010101 

2.630379373 

6 1 . 0010101801010101 

-2 . 60390166301 

62. 00101010101010101 

-1 . 335936042 

63 . 01010010101010101 

1 , 6.7'!39S39'^0i 

64 . 010101010010101VZ1 

“c 4 *4c! y-"Jc'c'yy 

65. 001001010000 

1 . 260073739 

66 . 01010101010100101 

1-1.0103151232 

67 . 00101010101010101 

-1.690710371 

63 . 010101010101010101 

0.477652776 

69 . 0101001010101010 

0 . ^07329263 

73 . 010101010101010101 

-1.075911357 

79 . 0100010101010101 

2 . 6'7 6353242 

30 . 0010101001010101 

1 . 03450 1 363 

3 1 . 00101010nZl01010 

-0 . 2 1 701745301 

32 . 00101001010101VZ1 

0.595362455 

S3 , 010101010101VZl0101 

-1 . 3296233301 

34 . 010101010101010101 

0. 3496392 19 

35. 0101010100100101 

-0 . 642044237 

36 . 01010101VZ101010101 

y 4 :■ t* 

37 . VZl010010101010101 

-0 . 720496367 

33 . 000010101010101 

-1 . 215329622 

39 . 0010101001010101 

-0 . 059234379 

93 . 01010101001010101 

0.441133365 

99. 0100010101010101 • 

“ 1 4 1 yT" 

001 . 00101010101010101 

2 . 00010106 1 701 


34951D-012 01 . " 0i4950i50l! 80i 0- Iyi000i0i01i 0i3 
DfiTR 



DToCP.FTE FCU'’'’’E'^ ^’D/\^fSFCP'V "C''JER SFECT^Ul'/' 
'-'OR UNEQUA-'^Y SEA^E TATA 


i-jyylt! ITlPL-lL-IT REnL^y !. H~H? 

060£0 C POICR SPECTRUM BY DISCRETE FOURIER TRhHSFORM METHOD 

00030 C NT = NUMBER OF TIME POINTS 

00040 C Ti;n ■= TIME POINTS.-. I = 1? NT 

00050 c :ki:i;i = i.'rlues of time series pit points ti:i:i...i = i, nt 

00060 C FST = STRRTING URLUE OF FREQUENOY 

00070 C FEND = ENDING URLUE OF FREQUENCY 

000S0 C FNIJ = NN = TOTRL NUMBER OF FREQUEICY POINTS IN SPECTRUM 

00090 C DELF = FREQUENCY INCREMENT 

00100 C FREQ = CRLCULRTED FREQUENCY 

00110 C SPFOUR = CRLCULRTED POI-iER SPECTRUM 

0Oil£0 DIMENSION T(500) 5 X(500) 

00130 PI = 3.14159367D0 

00140 RERD iSjli NT 

00150 1 FORMRTi;il0:i 

00160 FNT = NT 

0%1170 Du L I = 1 j NT 

001 130 ii! RERD US? S' .1 T 1 1 J ? Y f. I .( 

001'30 3 F0RMRTi:£F16.9:i 

00£00 HNORM = 0 

00210 DO 4 I = lit NT 

00220 4 '.^TIORM = XNORM + Xi:i:u«i«2 

00230 RERD i;3j5;i FST.* FEND.. FNM 

00240 5 FGR'MRTt:.3S15.S;i 

00250 DELF = (FEND - FST 3 .■••FNM 

00260 FREQ = FST 

00270 NhJ = FNI'J 

002S0 MRITE (4*63 FST* FEND* FNW* DELF* KNORM 

00290 6 FORMRT ( 5D 1 5 . 3 3 

00300 DO 91=1* NN 

(30310 SI = 01 . 

0101320 S2 = 0 

010330 DO 7 J = 1* NT 

00340 RRG = 2.0D0O-«PI*FREG!.'«T(.J3 

010350 CC = DCOS ( RRG 3 

00360 S3 = DSIN(RRG3 

003701 SI = SI + CC*X(.J3 

003S01 7 S2 = S2 -t- SS-^X(.J3 

0i0i3'?0i 'SPFOUR = 20101. 0iD0i* (SI + S2**23 (FNT*KNORM3 

010140101 I'jRITE (4?ii'.i FREL-!? sPFC'UR 

010410 S FORMRT (2F20. S3 

00420 9 FREQ = FREQ + DELF ; 

004301 9939 STOP 

00440 Ef-4T| 

END OF DRTR 
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■rfc- 


i KljT Fh'ijf'' 



5 1 I-0C 0 - 50495050B 

00 0- lylyyyO 

0 - 0104950150 

0.05091136 


0.03137588 

0.014S5149 

0.45515950 

01 . 0 1 9801 1 9S 

0 . 9;- 01013374 

0. 0c:4?5c4S 

3. 10084375 

0 . 0c:970iiE:97 

1 . 803:-i0439 

0^ 03465347 

01. 09096589 

01 . 0396039b 

1 . 80361300 

0 . 04455446 

0. 33318094 

01 ^ 0149501495 

1 . 31730030 

0. 05445545 

3. 13164406 

0 . 05940594 

0 . 50340859 

y j y t‘4-->--iS'4*4 

01. 0158601604 

01 ^ 0i6930i693 

y 4 1 ot* 

0.07435743 

1. 17118938 

0.07930793 

3.44393139 

0.03415S43 

9.56591945 

0. 039 10891 

15.95573054 

01 . 09401594 1 

8.98453135 

y - 0‘?Sy09'r'0 

3. 41890379 

01. 10139601401 

1 . 019463633 

0. 10891089 

0. 658901843 

0.11386139 

0.39891170 

0. 11881 188 

y< 5yl ;• 

0. 13376338 

0. 633079401 



y 4 i c!<:! 1 

1 . 08886bc4 

y 4 1 •*' '3 ^ ' l£' • ” 

6. 44894353 

0. 13861386 

14.84871381 

0. 14356436 

10. 50766744 

0. 14851485 


0. 15346535 

1 , 06355585 

0. 15841584 

0.33385471 

01. 16336634 

0.0110059801 

0. 16831683 

' 0.33561158 

0. 17336733 

3. 093580501 

0. 17831783 

3. 16817546 

0. 18316833 

0 . 46307088 

0. 18811881 

3.94337991 

0. 19306931 

0.-^5.336,50 

01. 198019801 

1.33765334 

0.30397030 

3.518678-4 

0.30793079 

01. 84306607 

0.3138:- 139 

0. 10678977 

0.31:- 831:- 8 

0. 03353864 


0.05195739 


0. 06'95498:- 

y - 

01 . 15416479 

0.33763376 

0.30465991 

0.3435-436 

0. 108:- 9933 

0.34753475 

0.87015305 

0.35347535 

01 . 3999 1 95 1 

0.35^43575 

01. 063;54450i 

01 . cbc:8 :- 6ii4 

01 . 013987:86.5 

01 . 86 8867 4 

0. 15733381 

y j 

0.37181631 

01 . 37733773 

01. 3:37551015 

0.38317833 

1.51631603 


48 


END 0 


0. E!S:‘ lc:S:-‘‘£ 

• - 4. l££95cQ2 


Pi . :-9E!0?9i=:i 

0„ 6?b931y4 

i'-’r '='• 

0, E9?0£9?l 

1,441 86665 


0.301980E0 

8, •44079083 


0 . :':06 -30 0 

0,97104549 


0.311SS119 

1,10936584 


0i. 31683109 

4,89183459 


0 j ijc! 1 yc! 1 >I« 

d 1 *!? 1 -,i'c=c!yo 


0 . 3c673c68 

01,49138389 


0,3316831? 

1 , 50443386 



0 , 4010:-; 1 4:-;6 


0,34158416 

8, 75189005 


0 , 34653466 

16,93447854 

. A' ■'* 

0,35148515 

80, 03493041 


0, 35643565 

4, 88869606 


0,36138614 

0,98387466 


01, 36633664 

14, 85985091 


0,3? 188? 13 

84 , 8467883y 


0,3?6c3?63 

6,70308788 


0,38118818 

5, 5345:33001 


0,38613863 

:3, 6:3708011:3 


0, 39108911 

1,18160418 

vV;’; "• 

0, 39603961 

01, 840157 666 


0 , 400990 1 0 

4,18341794 


0 , 405940601 

8, 84591403 


0, 41089109 

0, 14599885 

‘ 'If* : -'r 

0,41584159 

6, 43894989 


0,48079808 

81,118471:38 


0, 48574858 

8,86804158 ' 


0,4306930? 

0,80i45745i6 


0, 4356435? 

5,01868818 


0, 440594016 

0, 89754738 


0, 44554456 

y - 31 y 


0,4504'9505 

1 4 3 1 1 6yy 


0, 45544555 

8,938001988 


0,46039604 

1 , 87195889 


0 , 46534654 

01 , 84958559 


0, 47089703 

0, 08848718 


0,47584753 

01 , 8067:3755 


0,48019808 

01 , 11 58S6'0i6' 


0, 48514858 

8, 177:36784 


01, 490099011 

0,73315589 


0,49504951 

0, 19018:3616 

- 

01, 50010100010 

0, 496674 '94 


DRTfi 


. •• 
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COMPARISON PPOC^AM FC^ '^COKTNC AT C'IFFER~NT "C’'.'F° S~FCT-nM 


ES T T A T ES - ^'0 F F J EF ACK f-IA. M - TUK EY , 

p, D DT7T /i Ip Tf TTT?»r iPTf ^ 

FOU.^TFF T.FANSFCPM, AND ^EAST S3I.’A."E5 

(C^’.E "Y VA''"^CT< 


'T'U'5; Q'pHTo TV TTFc PAS-M'^ ^ 0 ) 


■00010 0. 

ULiiPAR I suh 

! PROGRAM 

FOR LOOKING AT 

*.-*ARIOUS POiCR SPECTRUM 

060801 C 

COMPUTATIO 

NAL METHi 

fsT'*- 


0!0i0;-;i.'i 

IMPLICIT R 

:EAL*8 (A- 

-H? 0-Z) 


06040 

DIMENSION 

T (20100 j ! 

■'■ 1 1. .1 ? i lEy !* t”' 


060*501 

DOUBLE PRECISION N 

1.200.1 j S i 200.! 


06i0i6'0i 

p T = Sft 1 4 1 

■5926' D01 



060i?0 C 

READ IN* TO 

TRL NUMB! 

:r of points nt 

Brill T I f‘ie SEP. I ES ?; 

06010:0 

READ i. 0: j i .1 

f'iT 



060901 

1 r URMAT 1 . 1 1 U 

) 



010 1 001 

pNJ = f'lT 




061 10 

DC: 2 I = i 

NT 



£61201 

2 KLAD 1 S' 5 S: 1 

T(I)? K 

( i ) 


06130 

3 FORMAT (2F1 

6.9'.i 




1.3.3140 C 
&11150 C 
03160 C 
03170 C 
061801 0 
06190 0 


06830 iJ 


ft I 


READ !M fi‘..iERAGE SPACING BETiCEN DATA POINTS BELT. STARTING 
FREQUENCY FSTj ENDING FREQUENCY FENIN TOTAL NUMBER OF FREQUEN 
POINTS FNH.> FRACTIONS OF THE flAYIMUM LAG NUMBER IN THE DATA 
I. NAYLAG.i FRACTIj FRACTSj FRAlTSj and FRAlT 4 FOR UOMRARISUN UF 
THE MODIFIED BLACKMAN- TUKEY SPECTRUM AND CORRELATIOr^ FUNCTION 
SPECTRUM I'JITH THE FOURIER? SPFOUR? AND TliO LEAST SQUARES ESTI 


062001 C 


MATES :3MELL0 (FERRAZ-MELLO) AND :3PLSQ i:i..!RNICEK) . 

06210 


READ (3? 4:1 BELT? F:ST? FEND? FNI'J 

062201 


READ ( 3 j 5 ) FR ACT 1 ? FRACT2 ? FRACT3 ? FRACT4 

062301 

4 

FORMAT (4D1U.8) 

06240 

C 

.J 

FORMAT (4F10,3) 

0612501 


NH = FNi'J 

06260 


NMl = .NT - 1, , 

062701 C 


REFERENCING THE TIME POINTS TO THE ORIGIN. 

062:30 


TM(i;i = O.ODO 

Kic‘90 

<■ 

DO 6 I = 1? NMl 

06:300 


DELA = T(I + i:i -Tin 

06:3101 

6 

TMd + lj = TMlIi + DELA 

06:320 C 


ROUNDING TIME POINTS TM(Ii TO THE NEAREST MULTIPLE 


Ai..'ERAGE TIME INTERNAL BELT 


06 : 3 ' 4 Gi 


DO 7 I = IjNT 

06:350 


hlTM = (TM(i:i.ftDELT;i + O. 5 D 0 

06 :s 60 i 

"7 

TM(i;i = MTM ■ 

00:3701 C 


OBTAINING THE NORM OF THE ! 



.■'dHuRM ~ 



iiij ft; 1 = 1 j NT 

06401 V 3 

zZ- 

ftNORM = ftNORM + Y(Ij »«2 

064101 


DO 9 I = 1 ;i 2001 

06420 


N(Ii = 01 . 01 D 01 

064:301 

Q 

; 3 Ri) = 01 . 0 D 01 

0614401 C 


CALCULATING SUMS OF LAGGED 

06450 


MAXLAG = TM(NT) - TMCl) + 

06460 


DO 11 I = 1 j NT 

06470 


DO 10 J = I j NT 

064:301 


ILAG = TMiJj - TM(i:i 1 . 0 ] 

0614901 


NdLAGi = NdLhG^ l.oDo 

065001 

i 

SdLAGi = SdLHGi + ftdi^^ft 

065 1 01 

1 1 

CONTINUE. 


1.0D0 


50 


ccfce 

•70 18- I r: 1 MPXLPG- 

0050 0 

TF I'M *' I ) . FO. 0 0’r0 ' GO TO 12 


N r T j = C; f T '! ..- f-J T '1 

:':rC=-0 

'c rOf.TIM!.E 

00?C0 

TELF = (FEr-4Ti _ ppT'i.-'t^NN 

.'0FTP 

F8T = DEIF 

00580 

FFEC = F8T 

■■ 0580 

PU 80 IFPFQ = MIJ 

:^080i" 

81 = 0 

00610 

88- = 0 

>■ 0680 

80: = 01 

00601' 

;^4 = fi 

|-,0g*0 

. 85 = 0 

C0 650 

86 = 01 

. 1 66i- 

87 = 0 


= J? 

067' 

8Mit = 01 

■06?0 

8M5 = 0 


XMERH = 01 

00674 

IiO 81 J = 1 j NT 

0067'" 

01 XMERN = 7MFRN + X(J) 

‘:0676 

7MERN = 7M<^PN- PFI.ORT(HT;i 

I'l' 68^0 

PO 10.. J = IrOT 

{•' 1,, t 'r*’-’’ 

F.J = J 

0:0.700 

•PPG = 8.. 0r0*P^*FPFQ*T(.J) 

,'•0.710 

CC = PC 08 '■ PPG '1 

10 780 

88 = P8IN(RRG) 

■*■0700 

81 =81 + CO**? 

0C7.-10 

88 = 88 + C0*88 

' >*"?50 

80' = 80 + 88**? 

il0?60 

P4 = pii + CO*'!'^ (...O' 

.:0770 

85 = 85 * S8*Xi J:i 

J-r ^ 

pf'14 = + cc'*(k(..i;i-7merh:i 

^.0770 

8M5 = 875 + 88i«(K(..!)-KMERM) 

■0780 

86 = 86 + CC 

00780 

57 = 57 + 55 


iesee 

C £:?0 

KCS-Jf 


0F900 
y IF 00 
.:1010 
'*' l 0 cF C 
i: 000 
01040 '*■ 
,-iP=e 


10 08 = 88 + X0!'i 
P0f'1£' = liFI, CRT I'NT) 

RIM? = 81 - (?4«W9..-R0f'1£:':i 
FRl_ Tl-' — 1 .••■'R0!f'1c' 

FRCTl = I S? - FpCT?'i.--RlMc- 
R8M? = 8? - i'??*»?--R0M£;i - FRCTl 


■-C860 

H01 = P8QFT ( 1 ,. 0iP0.-R0iM8;i 

10‘80'0 

R1 = D8QPT ( 1 . 0P01.--R1 M?) 

icsic 

PS = P8CPT ( 1 . 0ip0i..- R?f'1?;i 

:.0'rc0 

1^1 Cl = piW5r.i4 

rc900 

FRCT8 = 8? - ( (86*87). ••• 

IC'F-JC 

FRCTl = R1 *R2*C1 •*FPCT8 

•^0850 

C8 = R8*8M5 - FRCTl 

00860 

FRCT = C7**8 + C8**? 

:. 0'r70i C 

!■ Lh'Pbi; — l LU irF'lr C ! KUPi 

'■0080 

RMELLO = 1 00 . 0iP0i*FRCT.'' 


PET = 81-^83 - 88**8 
R = i 80*84 - 8?*85:i .xTCT 
.8 = i:-88*84 + Sl*85:t-'TiET 
i.iRMTCEK'V- LFPST SOURCE ESTIMRTE SPLSQ 
8PLSC = (100. 0r'0*i'84*R + S5*Bil ) .'-XMOPM 

ri8CFETE FCnpIEP TPRNSFOPM Fi0l-iEP SPECTPIJM E8TIMPTF CPFCI !p 
SFFOUP = 8vZie.0P0*( (84**8 + S5**8:i (FHT*yrJQPM:i 1 


OF LRG PPODUCTS EOURL TO FPRCTIOM8 QF MRXIMUM LRG 


1 060 

8ur'ii = 0-0P0 

1078 

8UM8 = 0C0P0 

• 080 

.5Ur.15 = 0,0PCi 

1 0'r 0 

5IJf1il = 0 ,, 0T|0 

1100 C 

setting Nl.if'lBI 
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01 1 1-C 


NLPGl 


FRACT 1 

pr 



^Jl2F 


NLAG2 

= 

FRACTS 

»MAXl. AC- 



1 1 'lip 


NLAG3 

= 

FRACT3 

:«MAXLAG 



oil 40 


NLAiB4 

= 

FRACT4 

^MAXLAij 



01150 


D*5 1 

C; T 

= 1? 

NLRG4 



01 160 


FI = 

T - 

1 




01170 


ARU = 

d 4 

0D0«PI 

4FRECf»F 

r«DELT 

01 1S0 


IFil- 

GT- 

NLAGl :* 

GO TO 

15 


011 •30 


SUNl 

= 

UMl + 

N *: I j «DC 

OS 

(arg; 

01200 

15 

if*:t- 

GT. 

NLAG2:* 

GO TO 

16 


01210 


SUM2 

= 

UN2 + 

S ( I j ■*DC 

OS; 

(ARC; 

01220 

16 

iF*;i- 

gtJ 

NLAGS:* 

*GO TO 

1 


012:50 


SUM3 

= 

UM3 + 

S f I ;* *DC 

OS 

(arg; 

01240 

17 

SIJM4 

= ■: 

IJM4 4 

S ( I 1 «DC 

OS 

(ARIj! 

01250 

IS 

CONTINUE 






ei£66 C MOrUFIEIi BLHCKMfiH-TUKEV SFECTF'JM I'JITH LRGS = "RFiCT 1 ^MFXLFiQ 
01E70 SFCGRl = E00. Glitters >0Ii0^SLlMl - .•■'r'riGRN 

01ES0 C CORRELATION METHOU SFECTRUM l-ilTH LAGS = FRACTE^MAXLAG 
01 £90 SFCORE = £G0. ODO^ (£. £iD0*SUM2 - XlCFNi ■■■■ iTNTSJTluRMJ 

01380 C CORRELATION METHOD SFECTRUM l-ilTH LAGS = FRACT3»MAXLAG 
01310 SFC0R3 = SOiOi.ODti^iSceiDOi^SUMS; - !-THORML-' (FNT*>TiORMj 

01320 C CORRELATION METHOD SFECTRUM I'nTH LAGS = FRACT4 «MA:;:LhG 
013315 SPC.I5R4 = 20OL 0Dti* (2. 8D0^SUM4 - :#CiRM FNT'^XNORM ) ’ 

013401 InR ITE 1.4? 19.1 FREQ? SFCuRl ? SF'CURL? SFL'URS ? SFC0R4 ? illELLO? SFFLiUR ? 

01350 ISFLSQ 

01360 19 FORMAT (F7., 4? SF9.. 4!i 

0*1 ’570 201 FRE*.;! = FRE*Ii + DELF 

01 ’330 9999 STOF 

01-390 END • 

END OF DATA 


ilCMFARISON OF NORMALIZED FOICR SFECTRUMS V-: FOICR) FOR THE MODIFIED BLACK- 
MAN-TUKEY METHOD? THE CORRELATION FUNCTION METHOD? FOURIER TRANSFORM AND 

LEAST SQUARES TECHNIQUES 

NOTE; THE SIMULATED DATA 'is'a TIME SERIES SAMPLED AT EQUAL INTERi..*ALS OF TIME 
(0,72 SECONDS APART) 0'..*ER A TIME SPAN OF 108 SECONDS- THE SERIES CAN BE 
. REFRESENTED BY THE FOLLONING EXPRESS ION- 

Xn;i = SIN(ARGl*Ii + i::OS(RRG£*I) + SIN (AR*33«I :* + COS *: ARG4« I i + SIN(ARG5*Ij 

l■JHERE AR.Gl = £-0»FIW0-37 
ARG2 = 2- 0»FI»0- 14 
AR*53 = 2-0**PI*«0-4£ 

AR|j 4 = £ - 0*P I *0 - 65 
ARI55 = 2-0»FI«l-09 

AND i'JHERE I kANGES* FROM 0* Tu '39-36 3E*-1JNDS- 'rT = 3- 141592fcF 


THE FIRST COLUMN IS FREQUENCY IN HERTZ 

THE SECOND iBOLUMN IS THE RELATI*.*E FOICR IN PERCENT? COMPUTED BY THE 
MODIFIED BLACKMAN-TUKEY FOICR SPECTRUM METHOD USING MEAN LAGGED PRODUCTS 
AND SUMMING UP TO 50X OF THE MAXIMUM LAG IN THE DATA- THE FOURIER 
TRANSFORM OF THIS MODIFIED AUTOCORRELATION FUNCTION IS THEN THE POiCR 
SFECTRUM - 

THE THIRD iJDLUMN IS RELATI*.C POICR IN PERCENT? COMPUTED E*Y ROUNDING THE 
TIME DIFFERENCES TO THE NEAREST MULTIPLE OF AN Ai..*ERAGE TIME INTERi.-*AL AND 
USING LAGGED OR CROSSED PRODUCTS UP TO 60X OF THE MAXIMUM LAG IN THE 
DATA- THEN THE FOURIER TRANSFORM OF THIS FUNCTION IS CALCULATED - 
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THE FOURTH COLUilh IS REL 
TIME DIFFERENCES TO THE 
‘JSSNU l—HUi_jED UK UFiiJSU’t-D 
DRTfl- Then ! he r-iJijRlLR 


^TIijE POl'iEF' 
CRREST MUL' 

■■Or-T‘i m*T.“ i :(“• 

. Uj_'Ui_- I L=r 


^RhU'b ijRM UF 


PERCENT? COMPUTED BY PCiUNDING 
LE of RN Ri'ERRGE TIME INTER'IRL 
^ OF THE MRYIMUM LRG IN THE 
I IS 'function is' CRilcULRTED- 



THE FIFTH COLUMN IS RELRTICE POICR IN PERCENT? COMPUTED BY ROUNDING THE 
TIME DIFFERENCES TO THE NERREST MULTIPLE OF RN Ri.CRRGE TIME INTERi.'RL RND 
USING LRGGED OR CROSSED PRODUCTS UP TO lOO'; OF THE MRKIMUM LRG IN THE 
DRTR. THEN THE FOURIER TRRNSFORM OF THIS FUNCTION IS CRLCULRTEDu 

THE SIXTH COLUMN IS RELflTI'.'E POICR IN PERCENT? COMPUTED BY THE METHOD 
OF S. FERRRZ-MELLG CALLED R ’DRTE OOMPENSRTE DISCRETE FOURIER TRANSFORM' 
i-JHICH i'jlLL B:E PUBLISHED SOOrL 

THE SE'JENTH COLUMN IS THE RELRTIi.C FOI’JEP IN PERCENT? COMFUTED BY THE FOURIER 
TRRNSFORM METHOD l-JITH THE OESER'.ED TIME POINTS. 


THE EIGHTH COLUMN IS THE RELRTIUE POICR IN PERCENT? COMPUTED BY RN OPTIMUM 
LERST SQUARES PROCEDURE (THE URNICEK METHOD I'JITHOUT 'CONSTITUENTS’? THRT IS? 
RUERRGES? L I NERR TRENDS? ETC. ) . 
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t-l 
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0459 
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0. 

0649 

0 . 0649 

0 

£i6*49 

0. 

46*96* 

0. 
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-0^4445. 

y - 

007 1 

0.0045 

0. 


0.0045 
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0. 

4746* 

-01. 

7451 

0. 8415 

0- 


0. £i360 

0. 

0360 

0 . 0360 

0 

0:360 

0. 

4796* 

-01. 

71 £6 

0. 0048 


0£01 

0. 0035 

0. 

0035 
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0 

00:35 

0. 

4846 

0. 
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Pi£49 

0.0817 

0. 

£i£17 

0 . 08 1 7 

0 
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0. 
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0- 
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0. 
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0. 0089 
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0S£1 
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-0- 
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0 
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0. 

50146 
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1345 
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0- 
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0.5190 

-0 . y637 

■ -0 . 4774 

6. 

0.0018 

0. 00118 

0.0018 

0.0018 

0.5346 

0.3843 

0.4707 
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Ci L- 1 

0.01 85 

0,01 F;Fi 

0.01 35 

y. yii96 

0.4957 



0 . 0015 


0,0015 

0.0015 

0 . -Iio 4b 

—0. 301.501 

—01 . 56 1 01 

0. 1036 

y - y 1 to 

0 . 01 68 

0,01 fc'S 

0.0168 

0.5396 

—01 . 50i96 

01 . 0i90i9 

—01 . 30183 

••1* . *i.'0' 1 

01. 010113 

0.001 3 

0.0013 

0.5446 

0 13654 
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0. 1691 

0.0317 

0. 0317 

0 . 013 1 7 
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0. 5496 

0.4568 

-0 . 3446 

0.0547 
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01 . 010 1 0 
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0 . 0379 
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—01 . 5589 

-01. 19001 

0 . 010018 

0 . 01008 

01 . 01001S 

01. 00018 

0.5745 

— 0 . 3'763* 

0.3931 
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0 . 015018 

0 . 05018 

01 . 0508 
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-0. 4336 

0.4931 . 
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01. 001017 

0 . 00107 

01 . 00017 
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0.5S45 

0.4434 
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01.0694 
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0.5945 

-01.5739 
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01 . 01974 
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01. 6045 
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0.4417 

0.7143 
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01 . 00018 
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01 . 54016 

0.3358 
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0 . 6395 
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7. 3163 
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0.6645 
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01.6761- 
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0 . 0470 
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0001G C DfiTfl GEHEPFlTOR FROGRRM 
000C0 IMPLICIT REFILLS i:R-HjO-Z:i 

0001301 l! I MENS I OM K C 5001 ) » T ( 500i J 

00040 n I MENS I ON • NT I M ( 500 ) ? RN 1 ( 300 ) 

01001501 D I MENS I ON TM ( 500 J 5 KM ( 500 ) j RN ( 500 ) 

010000 Ti I MENS I ON DEL (5001) 

00070 C SETTING UP SIMULRTED DRTR (T ( I ) j K( I j ) j I = 1 j NT i'iHERE 

0010S01 C Til) CONTRINS RRNDOM NOISE RNdi? ZERO MERN RND STRNDRRD 

000-?0 C DE'.'IRTION OF 0Lc'0i FROM R GRUSSIRN DISTRIBUTION 

00100 MT = 139 

00 110 FT = 3 . i 4 1 59367D0 

010130 RRGl = 3. 0D0 i-'*FT*0l37D0i 

00130 RRG3 = 3. 0D0*PI*0L 14D01 

00140 RRG3 = 3.0D0i*PI*0.43D0 

001501 RRG4 = 3^ 0D0i*Pr*0L65D0i 

0101601 RRG5 = 3.0D0*FT--«1.0i9D0i 

00170 DO 1 I = IjNT 

00 1 S01 T ( I !i = 01 . 73D0* 1 1 - 1 ) 

001 1 90 RN ( 1 3 = BRRN 1 ( - 1 j 0 ? 1 37S7 > 0 . 0D0 ? 0 - 30D0 ) 

0030101 T ( I j = T i: I j + RN ( I ) 

003 1 0 K i; I ;i = DS I N C RRG 1 «T f 1 3 3 + DCOS ( RRG3*T ( 13 3 + DS I N ( RRG3^«T i 13 3 + 

00330 1 DCOS { RRG4>T (13 3 +DS I N ( RRG5'-*T (13 3 

00330 STDE'J = 0L30D0*K(I3 
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00£4y 


RHl ( I i = ERRHl (-1 j 0? 13737? y« 0D0? STDEUj 

000501 C 


HDl'Ihij '-:RU55 IhN HuIoE I’ilTH rlEHN hNI 

tit10t« 0 


THE HMPLITUDE OF THE FUNCTIuri . 

00070 


H (13 = ;■< i I j + Rf -i 1(13 

00000 

1 

COrHTIHUE 

0100001 


riQ 3 I = 1?HT 

003010 

C. 

HTIM(I3 = 1 

01013101 C 


PUTTING "HOLES’ IN THE BfiTH 

001300 


liQ 334 I = 1?3 

0i0i:?30i 


NTIMi: 1+103 = 0 

01013401 


NT I M( I +40 3 = 0 

0101350 


NT I N i: 1+703 = 0 

0101360 

334 

NTIM( 1+303 = 0 

0103701 


L = 0 

00330 


DO 3 I = 1?NT 

0103301 


IF (NT IN ( 1 .1 .EQ. i;i L = L + 1 

00400 


IFi:NTIM(I3 .EQ.13 TN(L3 = 7(13 

010410 

3 

IF(NTIN(I3 .EQ, 13 XN(L3 = K(I3 

0104301 


NT = L 

0104301 


DO 51=1? NT 

0VJ440 

C 

NRITE (3? 63 TM(I3? :7M(I3 

0104501 

6 

FORMfiT(3F16.33 

010460 


STOP 

01014701 


END 




004301 


REhL FUNCT I ON EmRN 1 «3 ( I ? I KEY ? I FRN ? hMEhN 

004301 


I MFL I C I T REFiL-'^S ( R— H ? 0— Z ? $ 3 

00500 C 


3D.. -.THE DESIRED STRNDRRD DD.3 

003510 C 


RNERN...THE DESIRED MERN. 

005301 C 


H... THE POP SIZE. 

00530 


DRTfi IHERE.--137S7.-- 

010540 


DRTR H.-'36.0D0.-- 

01015501 


IF ( I KEY 3 5? 4? 4 

010156.01 

4 

I HERE = I FRN 

01015701 

5 

IF (13 6? 7? 7 

0101530 

6 

CRLL GROSS ( IHERE? SD? RNERN?'.3fiL? H3 

0105301 


I FRN = IHERE 

0101600 


GO TO 3 

0106 1 0 

y 

C RLL RRrHDU ( I HERE ? I FRN ? '.'RL 3 

006c0 


IHERE = I FRN 


0: 

BRRNl = URL 

01016401 


RETURN 

0106501 


END 

tiu660i 


SUERuUTINE gRUSS. ( IK? Sr RM? '.i?HJ 

001670 


I NPL I C I T F'ERLi*3 ( R-H ? 0-Z ? :T 3 



K = H 

MM63H 


R = 01 . 0D0 

00700 


DO 50 I = 1?K 

00710 


URLL .RntHDU 1. IX' IV' V i 

0y7c!0 


IK = lY 

00730 

501 

R = R + Y 

010740 


HO = H.-13. 0D0 

010750 


H3 = H--'3 . t.iD0i 

00760 


1.3 = (S«(R-H33 3 - DSQRKHCO + RM 



RETURN 

0017S0 


END 

001790 


SUEROUT I NE RRNDU ( I K ? I Y ? YFL 3 

0101000 


I MPL I C I T REPLIES ( R-H ? 0-Z ? 3 

001310 


DRTR JJJ5.--1037.-- 

0101330 


lY = IK«JJJ5 



IF ( IY3 5? 6? 6 

0013401 

c; 

lY = lY + 3147433647+1 


56 


8Ot850 

6 YFL = lY 



YFL — YFL*4 4656618D~':? 


00S70 

RETUF'N 



END 
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INPUT FOR CGMPRRISOH RUN WHERE THERE 
hRE holes RND £0:; RRNnOM NOISE ON TIME POINTS RNH Se 


10 ? 

4 £i?9 1 576 14 1 4 1 5946c£59 

0 4 99580 1484 1 4 8 1 8866086 


14388411039 

0.498414689 


—3 , ;-'0F;3636'5 1 

84 997585039 


34388048381 

■ 0 , 4'34'Jc;c;t.33 

44 546576938 

-1 4 709678913 

3 * cc!c'c! 1 1 c‘*43 

-1 4 463765533 

54641358839 

1 4 360980456 

64458689013 

0 4 859658 1 46 

13.019701019 

-0 4 86839 1 863 

134450376838 

-1 4 088413756 

14.466855535 

0.580768550 

15. 185540637 

04 796338959 

15.677851407 

—04 1S8540888 

1 6 4 46.6710930'^ 

—1 4 83449 101 

17. 548493564 

-04 '879947968 

18. 198036503 

-3.811410844 

1 S 4 600089 1 59 

-0.588001981 

19.380675883 

1 4 406833566 

80.077563841 

1.407787617 

80 4 93 r 83435 

-8.441510330 

81.690903398 

1.431447108 

88. 864556189 

3. 138538859 

83.477108488 

-8.670538876 

83.668130865 

-3.395948041 

84. 06I857784F1 

-04 1 16083536 

85. 886717966 

-1.918'834136 

86 4 £ 189 1 68837 

-0.897399918 

86.457930108 

-1.189143396 

87. 137409878 

0 4 033985866 

87.949730843 

1 4 533849944 

34.754894694 

-1 4 190770707 

35. 141514086 

1 4 688 1 1 4086 

36.085367381 

8.886819038 

36 4 900585064 

1.3915997V38 

37.648581910 

-8.559718441 

38 4 307 1 34884 

0.647067188 

38. 816997373 

1 4 706898984 

3R^ R0R438849 

— 1 , ;-;S4'^;-;i-:;30;-; 

39. 961538745 

-1 4 897087577 

40.915668131 

04 166577117 

41.943580713 

-1 4 384453588 

4c. 678473609 

-04 734808851 

43 4 503630064 

3.613903817 

43.714175443 

0 4 660006880 

•^■'4,734341150 

-1.701810978 


i:-; PNPLITUHE NO I 
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— b 4 ct'4c6'00 1 

-1 4 489468438 

46 4 8'860'8'?847' 

1 j 10 

4:^4 854866138 

-1 4 7773010554 

47 4 363837886 

-1 4 603746868 

4:-'4 :''fc6'8:’'466t^ 

04 594768387 

484 787141868 

1 4 8083081 1 

494 4701 '?9855 

1 4 883698818 

56.4 881 1'?8848 

-1 4 166704118 

56.4 651880109 

-1 4 198694774 

574846968851 

1 4 734877406 

58 4 8833346017 

01 4 7S6.66.84'?4 

594359984810 

-3 4 7010379578 

594683354489 

01 4 0185858369 

60 4 453'597687 

-0 4 46.'9S65 1 701 

601 4 7 68638478 

-04 786.8315101 

61 4 '?41 37 1890 

-04 ■?66816318 

684 40109'94344 

0 4 55016'? 77601 

634408638830 

14043751863 

634898181756 

-3 4 8'?867570iS 

704530775918 

04578765844 

7 1 4 86401401860! 

-04 0i81'?69185 

784 073'8575'8'? 

14014764318 

784 88447 10'37 

-0 4 704 189889 

40i0i5'r0S70i 

—01 4 S‘?3i60i5754 

744 81745.'5883 

-01 4 8014884.5018 

744880541358 

-84914847780 

754506811513 

-04 56445'? 188 

76.4 515097886 

84 4'?8'?99658 

76 4 94vli86vli487 

1 4 794'37'8737 

774 671875990 

-1 4 665761418 

78 4 7SiS358786 

84061889107 

784999363088 



79 . 650866773 0 * 3 1 7345843 

80 * 509300443 -3 4 887773 1 77 


y 1 ,« c!:'' 1 

04 393115'?46 


—8 4 3884483! 1 4 

884 6001 805'?3 

”c! 4 t' 1 :*'5y 

83 4 65 1 838056 

0 4 348034844 

884 S7'?94 1038 

04534444548 

844877768074 

1 4 01638018509 

854731179684 

01. 3'!!^5r36':!^620i 

86 4 4c'.:! 436 1 ti6 

84758509489 

864 978378891 

—01 4 886.6c!843o 


87 4 857 1 68940 - 1 4 808375389 

88 4 8 1 4035 74 -0 4 588388975 

89 4 593366947 - 1 4 1 1 0 1 48469 


'?04 141111110 

-04466807888 

90 4 877'?36809 

01 4 85'? 7870173 

9 1 4 564406 1 85 

-l4'?44701843 

'=!l4'^"1.5'Zi6.35,i3 

1 4 63605'? 130 

984 668 190844 

84 01787'?417 

^ rc!c' 1 1* : 

84056133478 

944816340853 

-1 4 '?67563774 


94 4 965 1 1 57 6c 4 44o6o6 114 

95 4 834780 1 99 0 4 884383534 


?6 4 84800101008 

“*c! .j oiz3yy5yuy 

?7 4 0164845663 

-34 148'959531 

?84 164480065 

84455408584 

?8 4 763961 196 

-04 189171878 

?9 4 543749434 

-14 0051060198 


04 01000 VZl4 :'0000 8004 0000 04 5000 04 6000 tl4 70!0£> 1 4 t!t100i '04 780.0 

END L'F DRT R 



COMPfiRISOH 
riHri— i iJKEV 


Uh r"'lJr!r‘lPL KUi'-iLR wPEL- ! FJjMo 

h-jrTHuIij THE CCRRELhT I OH FUhC^TIOM 
LEPST SQURRES “ECHN I 



-iE iluli I F‘ I EB BLF 




THE SIMULhTED DRTR IS R TIME SERIES l-HERE THE TIME POINTS NERE ORIGIMRL.LY 
SRMPLEB RT EQURL INTERURLS OF TIME ( 0.72 SECONDS fiPRRT ) 0'.iER R TIME SPRN 
OF 100 SECONDS. HOLES ICRE THEN PUT IN THE DRTR i. 32 OUT OF R ''OTRL OF 139 
POINTS MISSING OR REOUT 22R; J . R-RNDOM NOISE l-jlTH R GROSS I RN DISTRIBUTION 
HRUING ZERO MERN RND R STRNDRRD DE'.iIRTION OF 0.20 UNITS MERE PDDED TO THE 
REMfilNING TIME POINTS RND THE FUNCTION CRLCULRTED RT THERE PSEUDO RRNDOM 
POINTS. FINRLLVj RRNDOM NOISE FROM R GfiUSSIRN DISTRIBUTION HRlilNG ZERO MERN 
.RND STRNDRRD DEUIRTION EOURL TO 20F; OF THE RMPLITUDE OF THE FUriCTION i-iRS 
RDDED TO THE FUNCTIONRL i.jRLUES. (RRNDOM NOISE ON BOTH RXES!!. THE FREQUEr(CIES 
IN THE DRTR RRE 0.3? HZ? 0.14 HZs 0.42 HZ? 0.65 HZj RND 1.09 HZ. MITH THIS 
CHOICE OF FREQUENCIES RLIRSING IS EXPECTED RND CRN BE SEEN IN THE PRINTOUT 
BELOI-i. 


THE FIRST COLUMN IS FREQUENCY IN HERTZ. 

THE SECOND COLUMN IS THE RELRTIUE POMER IN PERCENT? COMPUTED BY THE 
MODIFIED ELRCKMRN-TUKEY POICR SPECTRUM METHOD USING MERN LRGGED PRODUCTS 
RND SUMMING UP TO 50X OF THE MRXIMUM LRG IN THE DRTR. THE FOURIER 
TRRNSFORM OF THIS MODIFIED flUTOCORRELRTION FUNCTION IS THEN THE F'OiCR 
SPECTRUM. 

THE THIRD COLUMN IS RELRTIi.iE POICR IN PERCENT? COMPUTED BY ROUNDING THE 
TIME DIFFERDICES TO THE NERREST MULTIPLE OF RN R 1 ..IERRGE TIME INTER 1 .JRL RND 
USING LRGGED OR CROSSED PRODUCTS UP TO SOX OF THE MRXIMUM LRG IN THE 

DRTR. THEN THE FOURIER TRRNSFORM OF THIS FUNCTION IS CRLCULRTED, 

THE FOURTH COLUMN IS RELRTI'.C POICR IN PERCENT? COMPUTED BY POUNDING THE 
TIME DIFFERENCES TO THE NERREST MULTIPLE OF RN Ri.JERRGE TIME INTERURL RND 
USING LRGGED OR CROSSED PRODUCTS UP TO 70X OF THE MRXIMUM LRG IN THE 

DRTR. THEN THE FOURIER TRRNSFORM OF THIS FUNCTION IS CRLCULRTED, 

THE FIFTH COLUMN IS RELRTIUE POICR IN PERCENT? COMPUTED BY ROUNDING THE 
TIME DIFFERENCES TO THE NERREST MULTIPLE OF RN RUERRGE TIME INTER'.iRL RND . 
USING LRGGED OR CROSSED PRODUCTS UP TO lOOX OF THE MRXIMHM LRG IN THE 
DRTR. THEN THE FOURIER TRRNSFORM OF THIS FUNCTION IS CRLCULRTED, 

THE SIXTH COLUMN IS RELRTI'.C POICR IN PERCENT? COMPUTED BY THE METHOD 
OF S. FERRRZ-MELLO CALLED R ’DRTE COMPENSATE DISCRETE FOURIER TRRNSFOF'M’ 
NHICA hJILL BE PUBLISHED SOON. 

THE SE'.CNTH COLUMN IS THE RELRTI'.C POICR IN PERCENT? COMPUTED BY THE FOURIER 
TRRNSFORM METHOD NITH THE OBSER'.CD TIME POINTS, 

THE EIGHTH COLUMN IS THE RELRTI'C POICR IN PERCENT? COMPUTED BY RN OPTIMUM 
LERST SAURRES PROCEDURE (THE '.ANICEK METHOD l■JITHOUT ’COr'STITUENTS’ ? THRT IS? 
R'.CRRGES? LINERR TRENDS? ETC.). 


0.0100 

-0.7933 

-0. 326S 

0.2249 

0. 0093 

0 . 0033 

0 , 0039 

0 , 0033 

0.0135 

-0 . 9097 

-0.0143 

-0. 5999 

0. 1933 

0.2147 

0, 1945 

0. 1932 

0.0169 

-0. IS'B'B 

0,7216 

0.0179 

0. 1332 

0. 1944 

0. 1736 

0, 1311 

0.0204 

0 , 900 1 

1 , 1 103 

1.4031 

0 - 609? 

0 , 5405 

0,5424 

0,5332 

0.023S 

1 , 5973 

1 . 0634 

1.7341 

1 . 7037 

1 . 4366 

1 , 6436 

1.4317 

0.0272 

1 . 5052 

1 , 0023 

0. 9306 

1 . 6553 

. 1.3756 

1 , 7427 

1 , 3334 

0. 0307 

0. 6304 

0,9931 

0. 4030 

0.4699 

0.5357 

0 , 5956 

0, 5362 

0.0341 

-0. 4433 

0,6300 

0. 5062 

0. 0062 

0.0131 

0, 0124 

0,0132 
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